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ABSTRACT 


An  initial  experimental  investigation  was  conducted  to 
examine  the  feasibility  of  NOj^  emission  control  using 
catalytic  reduction  techniques  in  the  jet  engine  test  cell 
environment.  A  modified  T-es  gas  turbine  combustor  and  an 
augmentor  tube,  21  feet  in  length  and  containing  a  perlite 
catalyst,  were  used  as  a  gas  generator  and  catalytic  reduction 
system.  Four  data  runs  were  made.  Three  runs  were  completed 
without  the  catalyst  installed.  Temperature  and  velocity 
profile  measurements  were  obtained  in  order  to  calculate 
augmentation  ratios  for  different  engine  fuel  to  air  ratios. 
NOjj,  CO,  and  unburned  hydrocarbon  concentrations  in  the 
exhaust  were  measured  to  provide  a  baseline  for  further  tests. 
A  fourth  data  run  was  made  with  the  perlite  catalyst  installed 
in  the  augmentor  tube.  A  64  percent  NOj^  reduction  was 
observed,  however,  the  large  pressure  drop  across  the 
catalytic  bed  deemed  the  current  configuration  impractical. 
Recommendations  for  alternative  configurations  are  presented. 
The  results  of  the  investigation  have  proven  that  further 
study  is  warranted. 
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I .  INTRODUCTION 


With  continually  increasing  worldwide  industrialization 
and  the  resulting  further  destruction  of  natural  cleansing 
agents,  environmental  concerns  and  pollution  control  are 
currently  of  paramount  importance  and  they  will  continue  to 
be  so  into  the  future.  A  major  portion  of  the  total  sources 
of  pollutants  are  the  particulate  and  gaseous  products  of 
combustion.  There  are  four  principle  chemical  classes  of 
these  pollutant  species  which  have  the  largest  impact: 

1.  Nitrogen  oxides  (NOj^) 

2.  Carbon  monoxide  (CO) 

3.  Organic  compounds  (unburned  or  partially  burned 
hydrocarbons  (HC)  or  soot) 

4.  Sulfur  oxides  (SO^^)  .  [Ref.  1] 

These  combustion  species  can  also  be  described  as  either 
primary  pollutants,  which  are  emitted  directly  to  the 
atmosphere  (soot  or  nitric  oxide,  NO)  or  secondary  pollutants 
which  are  formed  by  chemical/photochemical  reaction  of  primary 
pollutants  after  being  emitted  into  the  atmosphere  and  exposed 
to  sunlight  (i.e.,  ozone,  0^)  [Ref.  2].  Some  pollutants  fall 
into  both  categories  such  as  nitrous  oxide  (NO2) /  which  is 
emitted  directly  from  combustion  exhausts  but  is  also  formed 
in  the  atmosphere  photochemically  from  nitric  oxide  (NO) . 
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Of  the  total  nitrogen  oxides  (NO^^)  emitted  from  engine 
exhaust,  about  90%  is  NO  which  oxidizes  to  NO2  in  the 
atmosphere.  In  the  presence  of  sunlight,  NO2  decomposes 
photochemical ly  to  give  atomic  oxygen  (0)  which  combines  with 
molecular  oxygen  (O2)  to  give  ozone  (O2) .  The  reaction  of 
ozone  with  other  hydrocarbons  in  the  atmosphere  forms  carbon 
dioxide  (CO2)  and  other  secondary  products  which  result  in 
photochemical  smog  or  a  brownish  discoloration  of  the 
atmosphere.  [Ref.  3] 

The  U.S.  Navy  and  U.S.  Air  Force  have  been  very  interested 
in  studies  looking  for  economical  alternatives  in  controlling 
pollution  emissions,  specifically  soot  and  NO^^  emissions  from 
their  gas  turbine  engine  test  cells  (which  are  used  in 
research  and  engine  performance  recertification  after  a 
rework  or  overhaul) .  As  new  generation  engines  exhibiting 
somewhat  high  ar  exhaust  temperatures  and  pressure  ratios  are 
introduced,  it  can  be  expected  that  they  will  produce  more  NO^^ 
emissions  than  current  engines.  For  civilian  aviation,  such 
emission  standards  are  set  by  the  Environmental  Protection 
Agency  (EPA)  and  will  most  likely  become  more  stringent  in  the 
future.  The  overall  objectives  of  these  efforts  are  to  reduce 
the  emission  of  NO^^  and  other  pollutants,  minimize  the  damage 
to  the  atmosphere  and  satisfy  current  state  and  EPA 
regulations. 
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Commonly,  gas  turbine  engine  test  cells  are  comprised  of 
an  instrumented  static  test  stand,  an  augmentor  tube,  and  a 
vertical  exhaust  stack.  The  augmentor  tube  provides  an 
enclosure  to  reduce  the  velocity  and  temperature  of  the 
exhaust,  enable  potential  engine  noise  suppression,  and 
provides  a  capability  to  physically  "treat”  the  exhaust  gases, 
while  the  exhaust  stack  vents  the  gases  into  the  atmosphere. 

Currently,  emission  control  is  envisioned  to  use  some 
type  of  catalytic  reduction  process  within  the  augmentor  tube. 
This  could  involve  a  technique  combining  the  injection  of  some 
substance  (i.e.,  ammonia  or  isocyanic  acid  [Ref.  4,5,6])  into 
the  exhaust  flow  inside  the  augmentor  tube  and  a  catalytic 
surface  (i.e.,  iron  oxide,  perlite,  or  vermiculite  [Ref.  7]) 
installed  inside  the  tube  further  downstream.  This  technique 
of  NOjj  control  centers  around  the  chemical  treatment  of  the 
combustion  products  or  exhaust  of  the  engine.  The 
introduction  of  fuel  additives  in  the  engine  (thereby  changing 
the  combustion  species  in  the  exhaust)  and  the  use  of  staged 
combustion  have  also  been  studied,  but  will  not  be  addressed 
here. 

Previous  research  at  the  Naval  Postgraduate  School  since 
1982  has  primarily  focused  on  the  effects  of  fuel  composition 
and  additives  on  engine  exhaust  particulate  emissions  or  solid 
soot  concentrations  (as  last  reported  in  1988  by  Lindsay  [Ref. 
8]).  Lindsay's  research  utilized  an  actual  Allison  T-63-A-5A 
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engine  combustor  section  modified  to  allow  the  use  of 
nonintrusive,  optical  techniques  (three  wavelength 
transmittance/ forward  laser  light  scattering  measurements)  to 
measure  the  actual  size  and  concentration  of  solid  exhaust 
particles  (soot)  across  the  combustor  and  augmentor  tube. 

Conversely,  a  primary  goal  of  this  research  was  to 
investigate  an  effective  control  process,  entailing  a 
combination  of  exhaust  treatment  techniques,  to  accommodate 
the  operational  ranges  and  variables  experienced  in  the  gas 
turbine  engine  test  cell  environment.  The  challenge  exists 
in  determining  an  overall  control  strategy  since  catalytic 
reduction  has  been  found  to  be  very  temperature  dependent. 
In  addition,  the  technique  must  be  applied  without  severely 
compromising  the  proper  functioning  of  the  test  cell  (pressure 
drop,  augmentation  ratio,  etc.).  This  makes  it  particularly 
difficult  to  apply  to  different  engines,  engine  power 
settings,  and  corresponding  augmentation  ratios  in  the  gas 
turbine  engine  test  cell  environment.  Previous  research  at 
the  Naval  Postgraduate  School  has  not  touched  upon  this 
concept  before,  but  it  has  been  under  study  by  the  U.S.  Air 
Force,  at  the  Naval  Air  Propulsion  Center,  and  at  other 
government  laboratories  and  agencies. 

Recent  studies  have  observed  greater  amounts  of  CO  (700 
ppm  at  idle  to  45  ppm  at  75%  power)  and  unburned  hydrocarbons 
(950  ppm  at  idle  to  3  ppm  at  75%  power)  at  low  gas  turbine 
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power  settings.  At  high  power  settings,  larger  amounts  of 
(10  to  60  ppm)  and  smoke  have  been  observed.  The  general 
trend  reflects  that  high  concentrations  of  NOj^  are  not 
normally  present  together  with  high  levels  of  CO.  In 
afterburner,  gas  turbine  engines  may  emit  much  higher 
concentrations  of  both  CO  (900  ppm)  and  NO^^  (80  to  100  ppm) 
[Ref.  7]. 

The  variation  of  concentrations  with  power  settings  and 
other  variables  in  the  test  cell  environment  confirms  the 
difficulty  which  exists  in  providing  an  efficient  NO^^  control 
process  over  all  operating  conditions.  To  date,  studies 
indicate  that  NO^^  i  eductions  of  up  to  90-100%  are  possible 
depending  on  catalytic  bed  temperatures,  catalyst  volume  and 
composition,  and  engine  power  setting  [Ref.  7]. 

At  the  outset,  the  scope  and  following  objectives  were  set 
for  this  investigation  utilizing  the  Allison  T-63  combustor 
as  a  gas  generator. 

1.  Re-plumb  all  fuel  and  air  lines  and  re-instrument  the  T- 
63  engine  and  ensure  the  engine,  associated  hardware  and 
software  operate  satisfactorily. 

2.  Design  and  build  a  test  stand  for  the  21  foot  augmentor 
tube. 

3.  Install  and  calibrate  a  new  gas  sample  dilution  and 

conditioning  unit,  analyzer,  carbon  monoxide 

analyzer,  and  unburned  hydrocarbon  (UHC)  analyzer. 

4.  Determine  the  average  velocity,  temperature,  and  mass 
flow  rate  of  the  gas  exhaust  at  the  end  of  the  augmentor 
tube  prior  to  the  installation  of  the  catalytic  bed. 
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5.  Install  an  iron  oxide  (Fe-O^)  catalytic  bed  (available 
from  the  Von  Didier-Werke  dbrp. ,  West  Germany)  at  the  aft 
end  of  the  augmentor  tube  and  determine  the  effect  of  the 
catalytic  bed  on  NO^^,  CO,  and  UHC  concentrations, 
augmentation  ratio,  and  pressure  drop  as  the  combustor 
exhaust  passes  through  the  bed. 

6.  Measure  the  concentrations  of  NO^^,  CO,  and  UHC  before  and 
after  the  catalyst  while  varying* the  augmentation  ratio, 
engine  exhaust  temperature,  and  engine  fuel  to  air  ratio 
to  determine  the  practicality  of  using  such  a  technique 
in  the  gas  turbine  test  cell  environment. 
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II. 


EXPERIMENTAL  APPARATUS 


A.  COMBUSTOR 

A  full  scale  Allison  T-63-A-5A  Gas  Turbine  combustor 
(Figures  1  and  2),  as  modified  by  Grafton  [Ref.  9]  was  used 
to  generate  jet  exhaust  gases.  Grafton  [Ref.  9]  installed  a 
quench  manifold  forward  of  the  exhaust  nozzle  and  just  aft  of 
the  turbine  nozzle  block  to  simulate  the  temperature  drop  of 
the  combustor  exhaust  gases  which  would  normally  take  place 
upon  turbine  work  extraction.  In  this  investigation,  the 
quench  air  to  the  manifold  was  supplied  at  approximately  20 
deg.  F  (480  deg.  R)  and  at  a  flow  rate  of  between  0.5  and  0.6 
Ibm/sec.  The  required  quench  air  was  supplied  from  a  single 
air  line  branching  off  from  the  main  air  line  through  a  sonic 
choke  (Djjjp  =  0.2  37  in.)  sized  to  provide  a  0.5  Ibm/sec  flow 
rate  assuming  a  minimum  of  475  psi  air  pressure  was  provided 
upstream  of  the  choke.  The  quench  air  line  branch  was  located 
prior  to  the  main  air  sonic  choke  and  vitiated  air  heater. 
New  combustor  ignitor  and  air  heater  torch  transformers  were 
also  obtained  and  installed  prior  to  the  experiments. 

The  coitibustor  apparatus  used  in  this  experiment  was  moved 
from  the  test  cell  utilized  by  the  Lindsay  experiment  [Ref. 
8],  necessitating  installation  of  all  new  fuel  and  air  lines 
as  well  as  re-instrumentation  of  the  combustor.  Prior  to 
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initial  engine  test  and  data  runs,  all  fuel  and  air  lines  were 
leak  tested  under  high  pressure. 

Pertinent  T-63  engine  data  is  provided  below  in  Table  I 
at  standard  sea  level  static  conditions. 


Table  I.  T-63  PERFORMANCE  RATINGS  [Ref.  10] 


Rating 

f 

^t 

(deg  F) 

Takeoff 

0.019 

3.17 

0.061 

1380 

Military 

0.019 

3.17 

0.061 

1380 

Normal 

0.017 

3.04 

0.053 

1280 

90%  Normal 

0.017 

2.95 

0.049 

1226 

75%  Normal 

0.015 

2.82 

0.043 

1148 

Note!  compressor  ratio 
engine  length  = 
width  =  19.0  in, 

=6.25,  P^  =  92  psia 

40.4  in,  height  =  22.5  in, 
dry  weight  =  138.7  lb. 

B.  AIR  SUPPLY 

Compressed  air  for  the  combustor  and  quench  manifold  was 
provided  from  a  3000  psi  tank  storage  system  (Figure  3)  . 
Compressed  air  was  supplied  to  the  tank  system  using  two 
compressors  with  an  in-line  air  drier  system  to  remove 
moisture.  Air  flowed  from  the  tanks,  through  several  valves, 
to  a  dome  loaded  pressure  regulator  operated  from  the  system 
control  panel  inside  the  control  room.  The  dome  loaded 
pressure  regulator  provided  a  stable  pressure  to  the  main  air 
sonic  choke  (D_  =  0.42  in.),  which  when  instrumented  with  a 

Cl 
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Figure  3.  Schematic  of  Air  and  Fuel  Supply  Systems  (Adapted 
from  [Ref.  11]). 
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thermocouple  and  pressure  transducer  was  used  by  the  computer 
to  determine  the  air  mass  flow  rate  to  the  combustor.  The 
sonic  chokes  were  sized  to  obtain  an  approximate  1.7  -  1.9 
Ibm/sec  mass  flow  rate  of  primary  air,  0.5  Ibm/sec  of  quench 
air,  and  a  90-100  psi  combustor  chamber  pressure.  The  quench 
air  branch  was  located  after  the  dome  loaded  pressure 
regulator  and  prior  to  the  main  air  sonic  choke. 

C.  VITIATED  AIR  HEATER 

A  vitiated  air  heater  was  installed  between  the  main  inlet 
air  sonic  choke  and  the  inlet  of  the  combustor,  providing  a 
means  to  vary  the  combustor  inlet  temperature.  An  ethylene 
charged  ignitor  torch  was  used  to  start  the  hydrogen  fueled 
air  heater.  Make-up  oxygen  was  added  downstream  of  the  air 
heater  prior  to  entering  the  combustor  inlet.  This  oxygen 
was  added  to  replace  the  oxygen  burned  in  the  air  by  the 
vitiated  air  heater,  ensuring  thac  the  correct  molar/mass 
composition  of  inlet  oxygen  and  nitrogen  entered  the 
combustor.  The  torch  and  heater  gas  introduction  into  the 
main  air  and  vitiated  air  heater  was  controlled  via  solenoid 
valves  remotely  activated  from  the  control  room. 

Sonic  chokes  were  placed  in  the  lines  between  the  heater 
fuel  and  make-up  oxygen  gas  bottles  and  the  point  of  injection 
into  the  vitiated  air  heater.  The  chokes  were  sized  assum:ng 
an  air  neater  pressure  of  250  psi  and  hydrogen  fuel  to  air 
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of  0.0015,  yielding  an  air  heater  temperature  of  860  deg. 
Rankine.  Using  the  one  dimensional,  isentropic  flow 
expressions  with  fixed  properties,  the  heater  fuel  and  make¬ 
up  oxygen  sonic  choke  diameters  were  determined  to  be  0.035 
and  0.070  inches  respectively. 

D.  FUEL  SUPPLY 

A  remotely  controlled,  pressurized  20  gallon  tank  provided 
fuel  to  the  combustor.  Nitrogen  was  used  to  pressurize  the 
tank  using  a  remotely  controlled  dome  loaded  pressure 
regulator.  From  the  tank,  NAPC  #4  (Suntech  4)  fuel  was 
supplied  through  a  series  of  filters,  through  a  throttle  valve 
in  the  control  room,  into  a  turbine  flowmeter,  to  an  electric 
solenoid  shutoff  valve  and  into  the  combustor  (Figure  3)  . 
Fuel  flow  rate,  in  gallons  per  minute,  was  available  from  a 
digital  display  in  the  control  room  and  was  provided  to  the 
HP  data  acquisition  system  via  the  HP  3497  DACU  scanner 
channel  25.  The  displayed  flow  rate  was  the  result  of  the 
output  of  the  turbine  flowmeter. 

E.  THERMOCOUPLES  AND  PRESSURE  TRANSDUCER  INSTRUMENTATION 
Chromel-Alumel  (Type  K)  thermocouples  were  used  to  measure 

the  various  temperatures  in  the  high  pressure  air  lines  and 
combustor  as  well  as  inside  the  augmentor  tube.  The 
thermocouple  and  pressure  transducer  outputs  were  provided  to 
the  Hewlett-Packard  (HP)  3497  Data  Acquisition/Control  Unit 
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(DACU)  of  the  HP-3054A  data  acquisition  system  for  recording 
and  flow  rate  calculations.  The  following  provides  a  summary 
of  thermocouple  and  transducer  locations  and  their  associated 
data  acquisition  scanner  channel  number  used  to  input  the 
measurements  into  the  microcomputer: 

VARIABLE  HP-3497  DACU  CH  NUMBER 


P_  (main  air  pressure)  24 

cl 

P_  (combustor  chamber  pressure)  23 

Phf  (heater  fuel,  H2  pressure)  22 

Pho  (heater  make-up  oxygen  pressure)  21 

T_  (main  air  temp)  60 

Tcin  (combustor  inlet  air  temp)  61 


"^exl  (combustor  exhaust  upstream  of  quench)  62 
Tex2  (combustor  exhaust  downstream  quench)  63 
Tho  (heater  make-up  oxygen  temp)  64 
Thf  (heater  fuel  temp)  65 
Taugup  (augmentor  tube,  upstream  of  catalyst)  66 
^augdl  (augmentor  tube,  downstream  catalyst)  67 
'^augd2  (augmentor  tube,  downstream  catalyst)  68 


F.  AUGMENTOR  TUBE  AND  INSTRUMENTATION 

An  augmentor  tube,  21  feet  long  and  two  feet  in  diameter, 
(Figure  4)  was  used  during  the  experiment.  The  tube  was 
suspended  on  a  rigid  castor  system  mounted  on  a  heavy  stand 
constructed  of  four-inch  angle  iron.  The  castor  system 
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Figure  4.  Photographs  of  Augnientor  Tube  on  Stand 


enabled  the  augitientor  tube  to  be  slid  into  the  proper  position 
for  each  data  run  from  its  storage  position.  The  purpose  of 
the  txabe  was  to  capture  the  primary  jet  exhaust  of  the  T-63 
combustor,  mix  the  exhaust  with  inducted  ambient  air,  and 
enable  the  air  mixture  to  be  treated  by  the  catalyst  located 
ins.'’.:ie  the  augmentor  tube,  4.5  feet  from  the  exit.  Gas 
samples  re  taken  prior  to  and  after  the  catalytic  bed  to 
determine  the  effectiveness  of  the  catalyst  in  scrubbing  the 
exhaust  of  NO^^. 

A  variable  diameter  blocking  plate  was  attached  at  the 
forward  end  of  the  augmentor  tube  at  the  T-63  exhaust  (Figure 

5)  .  Varying  the  opening  diameter  of  the  tube  enabled  a 
variable  augmentation  ratio,  i.e.,  the  ratio  of  the  induced 
ambient  air  (air  drawn  into  the  tube  by  the  jet  exhausting 
into  the  tube)  flow  rate  to  the  T-63  exhaust  flow  rate  (Figure 

6)  .  This  in  turn  created  variable  amounts  of  colder  ambient 
air  induction  into  the  tube  to  viscously  mix  with  the  T-63 
exhaust  gases,  and  also  provided  the  required  cooler  catalyst 
inlet  and  augmentor  exhaust  temperatures.  For  this 
experiment,  a  five  inch  diameter  orifice  plate  v:as  attached 
to  the  blocking  plate. 

In  order  to  determine  the  total  mass  flow  through  the 
augmentor  tube,  an  average  velocity  and  temperature  near  the 
exit  of  the  tube  had  to  be  measured.  To  find  the  velocity, 
static  and  stagnation  pressures  were  obtained  30  inches  from 
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the  exit  plane  of  the  augmentor  tube.  Static  pressure  was 
measured  with  an  eighth-inch  static  port  installed  flush  with 
the  inside  of  the  tube.  The  stagnation  pressure  across  the 
tube  radius  was  obtained  with  a  modified  United  Sensor  KT-18- 
C/A-12-C  Kiel  probe  (Figure  7)  enclosed  within  a  0.375  inch 
stainless  steel  tube.  The  Kiel  probe  was  mounted  on  an 
electric  motor -driven  traversing  platform  (Figure  8)  mounted 
on  an  adjustable  stand.  The  Kiel  probe  was  introduced  into 
the  flow  through  an  orifice  in  the  wall  of  the  augmentor  tube. 
The  probe  could  extend  from  2.5  inches  from  the  wall  to  the 
tube  centerline  at  11.75  inches.  The  traversing  Kiel  probe 
was  controlled  by  a  DC  motor  speed  control  box  in  the  control 
room.  Low  pressure  transducers  on  the  adjustable  stand  were 
connected  through  Pacific  amplifiers  into  an  Omega  strip  chart 
recorder  for  recording  scatic  and  total  pressure  as  well  as 
probe  position  in  the  augmentor  tube. 

Figure  9  provides  a  schematic  of  the  augmentor  tube  and 
associated  instrumentation  including  thermocouple,  catalyst 
basket,  and  sample  gas  probe  locations. 

G.  CATALYTIC  BED 

A  stainless  steel,  circular  basket  was  manufactured  to 
enable  different  varieties  of  catalyst  species  to  be  inserted 
into  the  augmentor  tube.  The  enclosure  (Figure  10)  was  23.25 
inches  in  diameter  and  five  inches  wide,  producing  a  maximum 
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Figure  8.  Photographs  of  Kiel  Probe  on  Traversing  Stand 


Figure  9.  schematic  of  Augmentor  Tube  and  Instrumentation 
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of  1.2  cubic  feet  of  catalyst  bed  volume.  The  stainless  steel 
screen  mesh  containing  the  catalytic  material  was  made  of 
0.063  inch  wire  with  a  0.187  inch  mesh  spacing.  A  removable 
door  along  the  width  of  the  basket  provided  for  removal  or 
replacement  of  catalytic  material.  The  basket  was  machined 
so  that  it  could  be  slid  into  the  augmentor  tube  and  secured 
4.5  feet  from  the  tube  exit  plane.  Assuming  a  flow  velocity 
of  40  feet/sec,  minimum  residence  time  through  the  enclosure 
was  0.01  seconds. 

H.  GAS  SAMPLING  AND  ANALYSIS  EQUIPMENT 

Figure  11  shows  a  schematic  of  the  sample  gas  flow  path 
to  the  analyzers  after  being  collected  by  two  stainless,  1/4 
inch  sample  probes  in  the  augmentor  tube.  One  sample  probe 
was  positioned  upstream  and  one  downstream  of  the  catalyst 
enclosure.  A  three-way  solenoid  valve  connected  the  heated 
sample  line  from  the  downstream  sample  probe  in  the  augmentor 
tube,  the  upstream  sample  probe,  and  the  heated  sample  line 
to  the  gas  analyzers.  The  electric  solenoid  valve  was 
operated  from  the  control  room  via  a  toggle  switch  on  the  T- 
63  control  panel.  The  sample  lines  were  set  to  heat  the 
sample  to  a  temperature  of  275  deg.  F  to  prevent  water 
condensation. 
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Figure  11.  Schematic  of  Sample  Gas  Flow  Path 
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1.  Model  900  Heated  Sample  Gas  Dilution  and  Conditioning 

Unit 

Sample  „.;gas  was  drawn  at  about  1.3  SCFH  from  the 
stainless  sample  probes  in  the  augmentor  tube  through  a  heated 
sample  line  and  introduced  into  a  Thermal  Electron  Model  900 
Heated  Sample  Gas  Dilution  and  Conditioning  Unit  (Figure  12) . 
The  Model  900  blended  dry  dilution  air  at  a  20:1  dilution 
ratio  with  the  gas  sample  to  reduce  its  dewpoint,  preventing 
the  formation  of  condensate  at  room  temperature  [Ref.  13]. 
The  conditioned  sample  (output  of  the  Model  900)  was  then 
delivered  to  the  following  instruments  at  Seven  SCFH  and  ten 
psig  for  further  analysis: 

1.  Thermo  Electron  Model  lOAR  Chemiluminescent  NO/NOj^ 
Analyzer 

2.  Thermo  Electron  Model  48  Gas  Filter  Correlation  (GFC) 
Ambient  CO  Analyzer. 

2.  Model  lOAR  NO/NO^  Analyzer 

From  the  Model  900  Sample  Conditioning  Unit,  a  portion 
of  the  sample  gas  was  provided  via  teflon  tubing  into  the 
Thermo  Electron  Model  lOA  Rack-Mounted  Chemiluminescent  NO- 
Gas  Analyzer  (Figure  13)  for  continuous  measurement  of 
nitric  oxide  (NO)  and  nitrous  oxides  (NO+NO2  or  NOj^)  .  The 
Model  lOAR  was  capable  -of  measurement  ranges  from  2.5  to 
10,000  parts  per  million  (ppm)  and  had  a  sensitivity  of  0.1 
ppm. 
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Figure  12.  Photograpli  of  Model  900  Heated 
Scunple  Gas  Dilution  and  conditioning  Unit 
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The  basis  of  the  instrument  is  the  chemiluminescent 
reaction  of  NO  and  ozone  (O^)  or  NO+O^  — ►  NO2+O2  .  Light 
emission  results  when  excited  NO2  molecules  revert  to  their 
ground  state. 

To  measure  NO,  the  gas  sample  was  blended  with  ozone 
produced  by  an  internal  ozonator  in  a  reaction  chamber  within 
the  Model  lOAR  (Figure  14) .  The  resulting  chemiluminescence 
was  measured  through  an  optical  filter  by  a  sensitive 
photomultiplier  (PM) .  The  filter  and  PM  responded  to  light 
in  a  narrow  wavelength  band  unique  to  the  chemiluminescent 
reaction.  The  output  of  the  PM  was  linearly  proportional  to 
the  NO  concentration.  [Ref.  14] 

To  measure  NO^^,  the  sample  gas  was  diverted  through  an 
N02-to-N0  converter  in  which  the  nitrogen  dioxide  (NO2)  was 
thermally  converted  to  nitric  oxide  (NO)  for  subsequent 
measurement  via  the  chemiluminescent  process.  The 
chemiluminescent  response  in  the  reaction  chamber  to  the 
converter  output  mixture  was  linearly  proportional  to  the  NO^^ 
concentration  entering  the  converter.  [Ref.  14] 

3.  Model  810  Total  Hydrocarbon  Analyzer 

A  raw  portion  of  the  exhaust  gas  from  the  heated 
sample  line  was  provided  to  the  Thermo  Environmental 
Instruments  Model  810  Total  Hydrocarbon  Analyzer  (Figure  15) 
at  atmospheric  pressure  at  a  sampling  flow  rate  of  2500  ml/min 
(5.3  SCFH) .  The  sample  gas  for  the  Model  810  bypassed  the 
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Figure  is.  Photograph  of  Model  810  Total  Hydrocarbon 
Analyzer 
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Mode  1  900  due  to  the  high  dilution  ratio  of  the  Model  900  and 
the  low  concentrations  of  hydrocarbons  expected  in  the 
exhaust.  The  Model  810  measured  hydrocarbon  concentration 
ranges  of  0.1  to  10,000  ppm  with  an  accuracy  of  0.1  ppm 
utilizing  an  internal  hydrogen  flame  ionization  detector.  The 
detector  operated  by  ionizing  volatile  organic  compounds  using 
a  hydrogen  flame.  The  gas  sample  was  mixed  with  hydrogen 
prior  to  the  flame.  External  combustion  air  was  provided  for 
combustion  of  the  sample  in  the  flame.  Normally,  sample  flow 
and  hydrogen  flow  were  equal  and  combustion  air  was  five  to 
ten  times  greater.  Upon  combustion,  the  hydrogen  flame  burned 
the  organic  compounds  in  the  sample  to  generate  carbon  dioxide 
and  water.  Carbon  ions  were  also  formed  in  the  process.  This 
process  occurred  in  an  electrical  field  between  two 
electrodes,  one  near  the  hydrogen  flame  and  one  around  the 
flame  or  a  collector  electrode.  The  potential  difference 
between  the  two  electrodes  caused  movement  of  ions  to  one  or 
the  other  electrode.  A  small  ion  current  flowed,  was 
amplified,  and  provided  to  the  microprocessor  system  of  the 
Model  810.  [Ref.  15] 

The  Model  810  required  calibration  with  a  known  span 
gas  (known  concentration  of  hydrocarbons)  ,  zero  air  (less  than 
0.1  ppm  concentration  of  hydrocarbons),  combustion  air  or 
oxygen,  and  hydrogen. 
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4.  Model  48  GFC  Ambient  CO  Analyzer 

Another  small  portion  of  the  conditioned  sample  gas 
from  the  Model  900  Conditioning  Unit  was  drawn  off  and 
introduced  into  the  Thermo  Electron  Model  48  Gas  Filter 
Correlation  (GFC)  Carbon  Monoxide  (CO)  Ambient  Analyzer 
(Figures  16  and  17)  at  a  flow  rate  of  about  one  liter  per 
minute  at  atmospheric  pressure.  The  Model  48  was  capable  of 
measuring  CO  concentrations  from  0.1  to  1000  ppm  with  an 
accuracy  of  0.1  ppm  utilizing  non-dispersive  infrared 
absorption  techniques.  Since  infrared  absorption  is  a  non¬ 
linear  measurement  technique,  the  Model  48  transformed  the 
basic  analyzer  signal  into  a  linear  output.  This  was  done 
internally  by  the  Model  48  by  storing  the  calibration  curves 
in  computer  memory  and  subsequently  using  the  curves  to 
accurately  linearize  the  instrument  output  over  a  desired 
range.  An  internal  temperature  and  pressure  transducer 
provided  outputs  into  a  microcomputer  to  make  corrections  to 
instrument  output,  resulting  in  CO  concentration  measurements 
which  were  unaffected  by  changes  in  sample  gas  pressure  and 
temperature . 

The  sample  gas  CO  concentration  in  the  Model  48  was 
determined  using  an  internal  Gas  Filter  Correlation 
Spectrometer.  GFC  spectrometry  is  based  upon  the  comparison 
of  the  detailed  structure  of  the  infrared  absorption  spectrum 
of  the  measured  gas  to  that  of  other  gases  also  present  in  the 
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Figure  16.  Photograph  of  Model  48  GFC  CO  Analyzer 
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Figure  17.  Flow  Schematic  of  Hodel  48  CO  Analyzer 
[Ref.  16] 
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sample.  The  technique  is  implemented  by  using  an  input  of  a 
high  concentration  of  CO  (known  span  gas)  as  a  filter  for  the 
infrared  radiation  transmitted  through  the  analyzer  (GFC) . 
The  basic  components  of  the  GFC  CO  spectrometer  are  shown  in 
Figure  18.  Radiation  from  an  IR  source  is  chopped  and  then 
passed  through  a  gas  filter  alternating  between  CO  and 
nitrogen  due  to  rotation  of  the  filter  wheel.  The  radiation 
then  passes  through  a  narrow  bandpass  interference  filter  and 
enters  a  multiple  optical  pass  cell  where  absorption  by  the 
sample  gas  occurs.  The  IR  radiation  then  exits  the  sample 
cell  and  falls  on  the  IR  detector.  [Ref.  16] 

The  CO  gas  filter  acts  to  produce  a  reference  beam 
which  cannot  be  further  attenuated  by  CO  in  the  sample  cell. 
The  nitrogen  side  of  the  filter  wheel  is  transparent  to  the 
IR  radiation  and  produces  a  measure  beam  which  is  absorbed  by 
CO  in  the  cell.  The  chopped  detector  signal  is  modulated  by 
the  alternation  between  the  two  gas  filters  with  an  amplitude 
related  to  the  concentration  of  CO  in  the  sample  cell.  Other 
gases  do  not  cause  modulation  of  the  detector  signal  since 
they  absorb  the  reference  and  measure  beams  equally.  This 
means  that  the  GFC  system  responds  specifically  to  carbon 
monoxide.  The  sensitivity  of  the  Model  48  is  increased  to  0.1 
ppm  with  a  lower  detectable  limit  of  0.02  ppm  by  using 
multiple  pass  optics  in  the  sample  cell  leading  to  a  large 
path  length  or  improved  sensitivity.  [Ref.  16] 


36 


Figure  18.  Diagram  of  Model  48  6FC  Bpectrometer  [Ref.  16] 
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I.  CONTROL  ROOM 


A  main  control  room  and  associated  control  panels  were 
located  adjacent  to  the  test  cell  in  the  NFS  Combustion  Lab, 
Building  217,  NFS  Annex.  The  control  room  provided  a  secure 
and  quiet  space  to  control  and  visually  observe  the  experiment 
as  well  as  house  the  sample  gas  conditioning  unit,  gas 
analyzers,  and  data  acquisition  system.  Figure  19  provides 
a  photograph  of  the  layout  of  the  control  panel  including 
controls  for: 

1.  main  air  (red  guarded) 

2.  T-63  air  heater  torch  and  gases  (red  guarded) 

3.  T-63  combustor  ignitor 

4.  T-63  combustor  fuel  system  and  tank  vent 

5.  fuel  tank  pressure  and  hand  loader. 

Adjacent  to  the  panel  was  a  digital  readout  of  the  JF  fuel 
flowrate  in  gallons  per  minute  (GFM)  as  well  as  a  jet  exhaust 
temperature  readout  for  safety  backup.  The  main  air  pressure 
gauge  and  hand  loader  were  located  on  the  solid  fuel  ramjet 
(SFRJ)  air  heater  console  near  the  T-63  combustor  control 
panel . 

J.  DATA  ACQUISITION  AND  REDUCTION  SYSTEM 

A  Hewlett-Fackard  HF-3054A  automatic  data 
acquisition/control  system  located  in  the  control  room 
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Figure  19.  Photograph 


provided  test  data  acquisition,  reduction,  and  system  control 
fcr  the  experiment.  The  system  included  a  HP-3497A  data 
acquisition  and  control  unit  (DACU) ,  a  HP  3456A  digital 
voltmeter,  and  a  HP-9836S  microcomputer  with  a  HP-9153C  hard 
disk  drive  (Figure  20)  .  The  test  controlling  software  program 
(Appendix  A)  was  writter  in  HP  Basic  5.1  and  loaded  from  the 
HP-9153C  hard  disk  for  each  test  run.  The  software  was 
programed  for  transducer  calibration,  setup  of  gaseous  flows, 
data  acquisition  and  reduction  for  each  hot  data  run. 
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Figure  20.  Photograph  of  H-P  Data  Acquisition  and  Control 
System 
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Ill .  EXPERIMENTAL  PROCEDURE 

Prior  to  any  data  runs,  the  T-63  combustor,  augmentor 
tube,  traversing  probe  system,  gas  analyzers,  and  associated 
software  were  checked  and  tested  to  ensure  the  equipment 
operated  satisfactorily  after  installation. 

The  sample  gas  analyzers  were  calibrated  with  their 
respective  zero  and  span  gases.  In  accordance  with  Reference 
14,  the  Model  lOAR  analyzer  was  calibrated  with  a  zero  gas 
of  less  than  0.1  ppm  NO  and  a  span  gas  containing  220  ppm  of 
NO  in  nitrogen.  The  Model  48  CO  analyzer  was  calibrated  with 
a  zero  gas  of  less  than  0.1  ppm  CO  and  a  104  ppm  CO  span  gas. 
The  NOjj  and  CO  analyzer  calibration  gases  were  routed  through 
the  Model  900  Dilution  and  Conditioning  unit  prior  to  entering 
their  respective  analyzer  for  calibration,  since  the  actual 
sample  gas  would  also  be  conditioned  through  the  Model  900. 
Since  the  Model  900  diluted  any  sample  input  at  a  20:1  ratio, 
all  sample  gas  readings  from  the  Model  48  and  Model  lOAR 
during  a  data  run  were  multiplied  by  20  to  obtain  the  actual 
specie  concentration  in  ppm. 

The  Model  810  Total  Hydrocarbon  analyzer  was  calibrated 
with  a  zero  gas  of  less  than  0.1  ppm  of  methane  (CH^)  in  air 
and  a  span  gas  of  50  ppm  methane  in  air  in  accordance  with 
Reference  15.  Since  the  Model  810  automatically  initiated  a 
calibration  sequence  upon  start  and  at  prescribed  intervals, 
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the  zero  and  span  gases  were  permanently  plumbed  to  the  unit. 
Before  connecting  them  to  the  Model  810,  the  zero  and  span 
gas  flow  rates  were  set  externally  (to  0.53  SCFH)  using  two 
flow  meters  and  needle  valves.  Due  to  the  calibration 
requirements,  the  T-63  exhaust  sample  gas  was  introduced 
directly  into  the  Model  810,  bypassing  the  Model  900.  All 
analyzers  and  the  Model  900  were  energized  one  hour  prior  to 
a  hot  data  run  to  ensure  proper  warmup  and  operation.  This 
normally  allowed  enough  time  for  the  Model  900  chamber 
temperature  to  reach  180  degrees  F  as  required. 

All  pressure  transducers  required  calibration  prior  to 
data  collection.  The  computer  program  *'T63NOX''  (Appendix  A) 
was  used  to  complete  the  calibration.  New  calibration 
constants  and  zero  values  were  obtained  and  entered  into  the 
program  for  data  reduction. 

The  Kiel  total  pressure  probe  and  static  pressure  port 
transducers  were  calibrated  with  a  manometer  and  strip  chart 
recorder  in  the  control  room.  The  calibration  resulted  in  a 
constant  of  10  mV  per  0.58  cm  of  water  for  the  static  port  and 
10  mV  per  0.53  cm  of  water  for  the  Kiel  total  pressure  probe. 

With  the  probes,  transducers,  and  analyzers  set,  the  Run 
Checklist  (Appendix  B)  was  completed  and  the  various  flow 
rates  (main  air,  air  heater  fuel,  heater  oxygen)  were  set  in 
accordance  with  the  ''T63NOX''  computer  program  on  the  H-P 
microcomputer.  On  cue  from  the  computer,  the  air  heater  and 
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combustor  were  ignited  from  the  control  panel  in  the  control 
room.  The  combustor  was  operated  at  normal  engine  operating 
fuel  to  air  ratios  (f  =  0.017-0.019).  The  Kiel  probe 
traversed  across  the  augmentor  tube  exit  plane  and  gas 
samples  were  taken  from  both  collection  probes  in  the 
augmentor  tube  and  processed  by  the  gas  analyzers. 

Originally,  four  hot  data  runs  were  to  be  completed.  Two 
runs  each  were  to  be  made  both  with  and  without  a  catalytic 
bed  inserted  in  the  augmentor  tube;  one  of  these  runs  without 
the  vitiated  air  heater  activated.  The  runs  without  the 
catalyst  were  made  to  establish  the  base  concentration  of  NO^^, 
CO,  and  UHC  and  to  determine  the  augmentation  ratios  from  the 
average  temperatures  and  velocities  across  the  tube  exit 
plane.  CO  and  UHC  levels  were  measured  since  it  had  been 
shown  (Ref.  4)  that  their  presence  can  significantly  change 
the  effectiveness  of  a  catalyst  for  NO^^  reduction.  For 
example,  NO  can  be  reduced  by  CO  using  a  rhodium  catalyst  to 
produce  CO2  and  N2 .  The  other  two  data  runs  were  to  be  made 
with  an  iron  oxide  catalyst  inserted  in  the  augmentor  tube 
following  the  work  of  Reference  6.  When  the  iron  oxide 
catalyst  became  unavailable  (due  to  the  supplier  deciding  that 
it  would  not  work  adequately  in  the  augmentor  tube 
environment) ,  a  vermiculite  catalyst  (an  aluminum-iron- 
magnesium  silicate)  was  ordered  to  evaluate  the  recent  results 
of  Reference  7  within  the  augmentor  tube  environment.  This 
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material  was  not  delivered  in  time  for  evaluation.  Course 
perlite  (an  inexpensive,  amorphous,  sodium-potassium-aluminum 
silicate)  was  then  obtained  for  testing  (also  based  on  the 
recent  results  of  Reference  7) ,  but  the  grain  size  proved  to 
be  too  small  for  the  catalytic  bed  screen  mesh.  A  finer  mesh 
screen  was  therefore  attached  to  both  sides  of  the  coarse  wire 
basket.  A  fourth  data  run  was  then  completed  with  the 
catalyst,  but  without  activating  the  vitiated  air  heater.  The 
run  was  made  to  evaluate  the  effectiveness  of  the  perlite 
catalyst  for  reduction  and  to  determine  the  effect  of  the 
catalytic  bed  on  the  augmentation  ratio. 
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IV.  RESULTS  AND  DISCUSSION 

A  summary  of  the  results  of  the  four  data  runs  is 
presented  below  in  Table  II.  A  more  comprehensive  data  output 
for  each  run  appears  in  Appendix  C. 

Table  II.  SUMMARY  OF  RESULTS 


MEASUREMENT  RUN  NUMBER 


1 

2 

3 

4 

air  heater  used 

no 

yes 

no 

no 

catalyst  installed 

no 

no 

no 

yes 

Tsi^  ,  (Ibm/sec) 

1.92 

1.95 

1.83 

1.91 

f 

0.017 

0.009 

0.019 

0.018 

Vavg  (ft/sec) 

41.5 

33.0 

36.2 

* 

AR 

1.17 

0.69 

0.78 

negl . 

Taugup 

903 

899 

997 

989+ 

^exl 

1550 

1374 

1693 

1555+ 

NOjj  (ppm) 

26 

14 

30 

50/18 

CO  (ppm) 

360 

82 

220 

520/220 

UHC  (ppm) 

135 

250 

58 

38/18 

Notes:  For  Run  #4: 

*  unable  to  calculate 
+  data  taken  early  in  test 
gas  concentrations  : 
upstrearo/downstream  of  catalytic  bed 
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It  is  evident  from  the  data  of  runs  1,  2,  and  3  that  as 
the  combustor  temperature  increased,  exhaust  levels  also 
increased  and  UHC  concentrations  decreased  as  expected  for  gas 
turbine  combustors  and  confirming  the  trend  found  in  Reference 
7.  CO  levels  were  quite  high  for  the  first  three  runs  except 
for  run  #2,  which  is  suspect.  As  the  combustion  temperature 
and  augmentation  ratio  decreased,  CO  levels  should  have 
increased.  This  behavior  was  not  observed  and  could  have  been 
the  result  of  either  the  low  combustor  flow  rate  or  improper 
functioning  of  the  CO  analyzer.  Overall  CO  and  UHC  levels 
were  somewhat  higher  than  reported  in  Reference  7 . 
Differences  could  have  been  attributed  to  defined  power 
settings  and  combustion  temperatures  in  the  engine  data  of 
Reference  7.  Additionally,  the  Model  810  UHC  analyzer  used 
in  the  tests  had  experienced  calibration  problems  during  the 
experiments. 

The  average  velocities  near  the  augmentor  tube  exit  (18.5 
feet  downstream  of  combustor  exhaust)  for  the  first  three  runs 
were  fairly  close  in  value.  The  velocity  and  temperature 
profiles  were  fairly  uniform  with  a  relatively  constant 
velocity  (constant  stagnation  pressure)  in  the  center  of  the 
flow  as  expected.  The  flow  did  exhibit  slightly  greater 
velocities  near  the  wall  of  the  tube,  possibly  the  result  of 
the  mixing  process  between  the  engine  exhaust  and  augmentor 
air  in  the  large  diameter  augmentor  tube.  The  lower  velocity 
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and  augmentation  ratio  during  run  #2  (vitiated  air  heater 
activated)  was  a  result  of  an  initial  combustor 
overtemperature  condition  at  the  normal  operating  fuel  to  air 
ratio  upon  ignition.  This  required  a  reduced  fuel-to-air 
ratio  to  be  used.  The  resulting  correction  was  too  severe  and 
produced  a  much  lower  combustor  exhaust  temperature. 

Run  #4  was  accomplished  with  a  perlite  catalytic  bed 
inserted  4.5  feet  from  the  txabe  exit.  Temperatures  entering 
the  bed  were  initially  989  deg.  R  (529  deg.  F) .  Due  to  the 
large  pressure  drop  caused  by  the  catalytic  bed  obstructing 
the  flow  through  the  tube,  the  augmentation  ratio  was  greatly 
reduced,  resulting  in  the  tube  temperature  being  increased 
substantially  as  the  run  progressed.  The  velocity  profile 
downstream  of  the  catalyst  exhibited  areas  of  flow  reversal 
and  substantially  higher  velocities  near  the  wall  due  to  an 
inexact  seal  between  the  wall  and  catalyst  enclosure. 
Additionally,  fine  particles  of  perlite  were  inducted  into  the 
Kiel  probe  head,  hampering  efforts  to  obtain  a  reliable 
average  exit  velocity.  Based  upon  the  measured  data  and  the 
observed  very  high  augmentor  tube  wall  temperature,  the 
augmentation  ratio  was  known  to  be  very  small. 

The  perlite  catalytic  bed  did  provide  for  a  64%  reduction 
in  NOj^  which  was  on  the  order  found  in  Reference  7  (for  a 
subscale  test)  for  a  perlite  +  MgSO^  bed  at  860  deg.  R.  The 
long  residence  times  (probably  on  the  order  of  0.01-0.03 
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seconds) ,  large  exposed  surface  area  of  the  catalyst 
particles,  and  higher  bed  temperatures  probably  contributed 
to  the  apparent  success  of  the  perlite  catalyst.  The  catalyst 
also  appeared  to  cause  a  decrease  in  CO  and  UHC.  This  may  be 
the  result  of  the  oxidation  of  CO  over  the  catalyst. 

It  was  apparent  from  these  initial  tests  that  the 
catalytic  reduction  of  NO^^  in  test  augmentor  tubes  is 
practical.  However,  to  not  adversely  affect  the  augmentor 
flow  rate,  the  catalysts  will  have  to  be  distributed  in  such 
a  manner  to  reduce  flow  rate  resistance.  Two  such  methods 
are;  (1)  placing  the  catalyst  material  in  tubular  structures 
(honeycombs)  and,  (2)  treating  only  the  central  region  of  the 
flow  where  the  NO^^  levels  are  at  their  highest  values. 

In  this  initial  investigation,  the  augmentor  tube  diameter 
was  made  large  and  the  augmentor  inlet  orifice  small  in  order 
to  provide  low  velocities  (30-40  ft/sec)  through  the  catalytic 
bed.  This  was  done  to  determine  the  effectiveness  of  the 
catalyst  under  the  most  ideal  conditions  and  with  a  minimum 
of  catalytic  material.  Augmentation  ratios  as  high  as  five 
and  velocities  as  high  as  1000  ft/sec  are  not  uncommon  in 
full-scale  test  cells.  Once  the  most  effective  catalyst  is 
found,  it  must  be  evaluated  over  the  full  range  of  velocities 
and  temperatures  encountered  in  the  test  cell  environment. 

Another  issue  which  must  be  addressed  is  the  impact  of 
soot  (which  can  accumulate  with  run  time)  on  the  effectiveness 
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of  the  catalyst. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


Most  of  the  initial  objectives  for  this  investigation  were 
met.  The  T-63  combustor,  associated  instrvimentation,  and 
software  operated  satisfactorily.  A  test  stand  for  the 
augmentor  tube  was  designed  and  built.  The  sample  gas 
analyzers  were  set  up  and  calibrated  and  a  sample  flow  system 
was  constructed  for  transporting  the  exhaust  sample  from  the 
probes  in  the  augmentor  tube  to  the  analyzers.  A  long 
response  time  was  noticed  during  the  tests  due  to  the  long 
sample  hose  lengths  required  to  transport  the  sample  to  the 
analyzers  and  small  orifice  size  in  the  three-way  solenoid 
valve  at  the  augmentor  tube  stand.  Movement  of  the  analyzers 
to  the  test  cell,  increasing  the  orifice  of  the  solenoid 
valve,  and  increasing  the  sample  probe  size  to  3/8  inch 
stainless  steel  could  minimize  sample  travel  time  and  distance 
to  the  analyzers  and  suction  resistance  in  the  tubing. 

The  general  velocity  and  temperature  profiles  across  the 
augmentor  tube  were  determined  to  be  relatively  constant  and 
used  to  obtain  augmentation  ratios  when  practical.  The 
perlite  catalytic  bed  resulted  in  a  64%  removal.  This 
supports  a  conclusion  that  perlite  is  a  viable,  inexpensive 
catalyst  material  which  could  be  used  in  the  test  cell 
environment.  The  large  pressure  drop  observed  across  the 
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catalytic  bed,  however,  makes  the  current  configuration 
impractical.  Alternative  configurations  might  include: 


1.  For  total  exhaust  flow  treatment,  a  honeycombed  catalyst 
enclosure  could  be  manufactured  which  would  allow  a 
larger  flow  velocity  through  the  bed  and  lower  pressure 
drop.  Systems  of  this  type  are  discussed  in  Reference 

6. 

2 .  Construction  of  an  extractive  system  involving  placement 
of  a  smaller,  coarse  catalytic  bed  centered  along  the 
augmentor  tube  center  axis  where  the  greater 
concentration  of  might  be  expected.  To  construct 
this  configuration,*  velocity,  temperature,  and 
concentration  profiles  would  have  to  be  found  along  the 
augmentor  tube  to  optimize  exact  catalytic  bed  placement 
and  NOj^  removal. 

Since  time  constraints  precluded  a  more  in-depth 
investigation,  future  efforts,  utilizing  the  current  apparatus 
and  various  species  of  catalytic  material  (Fe202  and 
vermiculite  when  available)  ,  should  concentrate  on  varying  the 
augmentation  ratio  and  fuel  to  air  ratio  and  measuring  NO^^, 
CO,  and  UHC  levels  across  the  catalyst.  The  effects  of  these 
variables  in  consonance  with  ammonia  or  cyanuric  acid  addition 
techniques  (as  discussed  in  References  4,  5,  and  6)  on  the 
effectiveness  of  the  catalytic  process  should  be  investigated. 

In  addition,  the  effects  of  sooting  and  much  higher 
augmentor  tube  velocities  and  augmentation  ratios  must  be 
evaluated  before  a  practical  solution  can  be  attained. 
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APPENDIX  A 


HP  BASIC  "TeSNOX"  COMPUTER  PROGRAM 


30 


BO 

90 

100 

US 

120 

130 

MO 

158 

160 

170 

180 

190 

200 

210 

Hi 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

a 

470 

460 

490 

500 

510 

520 

530 


9,  JAN  1990 

1  T63  COMBUSTOR  DATA  ACQUISITION  AND  REDUCTION  PROGRAM 

I  BMamE 

i  III  mSs WL"fw«iio« 

•  (4)  THE  TEST  SEQUENCE  AND  DATA  COLLECTION 

cvwoni  definition 

ANALOG  CHANNEL  NUMBER 
THROAT  AREA, AIR  FLO«  SOfTlC  CHOKE,  SQ  in 
THROAT  AREA, HEATER  FUEL  SONIC  pDKE,  SQ. IN. 
THROAT  AREA ’heater  OXYGEN  SONIC  CHOKE,  SQ.IN. 
BYPASS  AIR  FLOURATE 
DISCHARGE  COEFFICIENT,  AIR  SONIC  CHOKE 
DISCHARGE  COEFFICIENT,  HEATER  FUEL  SONIC  C^KE 
DISCHARGE  COEFFICIENT.  HEATER  02  SONIC  CHOKE 
AIR  SONIC  CHOKE  DIAMETER 
BYPASS  AIR  SONIC  CHOKE  DIAMETER 

AIR^HEATER  FUEL^^ONIC  CHOKE  DIAMETER 
air  HEATER  OXYGEN  SONIC  CHOKE  DIAMETER 
FUEL  IDENTIFICATION 
32.174 


!  SYMBOL 
!  A 
! 
i 
i 
i 
i 
i 
i 
i 

j 

i 

i 
1 

j 
i 


Aafr 
Ahf 
Aho 
BFair 
Cdair 
Cdhf 
Cdho 

Dairchoke 
Dbpchoke 
DateO 
Dhfchoke 
Dhochoke 
F«e]id$ 

.  Gc 
Gc=32,174 
!  Heaterftiel 

I 

\  Knair 
!  Kfifuel 
!  Knhf 
!  Knho 

!  KP.. 

!  Hair 
!  Maird 
!  Hfuel 
!  Mfueld 
!  Mhf 
!  Mhfd 
!  Mho 
!  Hhod 
a 

bpa 
Pbar 
Pc 
Phf 
Pho 
Ta 
Tbpa 
Tauaup 


heater  fuel  IDENTIFICATION 

AIR  SONIC  CHOKE  FLOW  RATE  CONSTANT 
FUEL  FLOW  METER  RATE  CONSTANT  (GPM/VDLT 
HEATER  FUEL  SONIC  CHOKE  FLOH  RATE  CONSTANT 
HEATER  02  SONIC  CHOKE  FLOR  RATE  CONSTANT 
PRESSURE  TRANSDUCERS  CONSTANT  (PSI/VOLT) 

AIR  FLOQ  RATE,  LBH/SEC 
DESIRED  AIR  FLON  RATE,  LBM/SEC 
FUEL  FLOy  RATE,  GPU 
DESIRED  FUEL  FLOW  RATE,  GPM  ^ 

HEATER  FUEL  FLOW  RATE,  LBH/SEC 
DESIRED  HEATER  FUEL  FLOW  RATE  LBM/SEC 

f«:  Sk\,  psia 

BAROMETRic  PRESSURE,  PSIA 
PRESSURE,  COMBUSTION  CHAMBER,  PSIA 
PRESSURE  HEATER  FUEL  SONIC  CHOKE.  PSIA 
PRESSURE  HEATER  OXYGEN  SONIC  CHOKE,  PSIA 
TEMPERATURE,  AIK  SONIC  CHOKE,  R 

IISI5SI;  SMffSpSB'cATALrsi. . 
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!  Taugdl  '  TEMPERATURE r  AUGMENTOR  TUBE  DOWNSTREAM  CATALYST,  R 

!  Testno*  TEST  I, T  NO.  ' 

!  Thf  TEMPERATURE,  HEATER  FUEL  SONIC  CHOKE,  R 

!  Tho  TEMPERATURE,  HEATER  02  SONIC  CHOKE,  R 

!  Tcin  TEMPERATURE,  COMBUSTOR  AIR  INLET  (HEATER  OUTLET),  R 

!  Tcind  TEMPERATURE,  DESIRED  COMBUSTOR  AIR  INLET,  R 

!  Text  TEMPERATURE,  COMBUSTOR  EXHAUST  UPSTREAM  OF  OUENCH,  R 

'  Tex2  TEMPERATURE,  COMBUSTOR  EXHAUST  DOWNSTREAM  (F  QUENm,^ 

BEEP  1000,.! 

PRINT  USING  'g' 

PRINT  USING  ‘br 
PRINT  •  T63  DATA  ACOUISITION  ' 

PRINT  USING  ’S/’ 

PRINT  'TURN  THE  PRINTER  ON  LINE* 

CLEAR  709 
CLEAR  722 

!  THE  RECORDED  VARIABLES  (VOLTAGES)  AND  LOCATIONS  ARE: 

!  (NOTE:  THE  MAXIMUM  ALLOWABLE  VOLTAGE  INTO  THE  SYSTEM  IS  1.2  VOLTS) 


I  VARIABLE  3A97  DACU  SCANNER  NUMBER  0__ 

Fa -  - 

'  Pbpa  24 

Pc  23 

Mfuel  25 

Phf  22 

I  Pho  21 

Ta  60 

Tcin  (inlet  air)  fel 

Te»l  (upstrean  of  quench)  62 

Tpx2  (dounstrean  it  qgench)  63 

Tho  64 

Thf  65 

Taugon  66 

Tauqdl  67 

Taugd2  6B 

ALL  Plow  rates  are  calculated  using  the  ONE-Dln£i, signal,  ISENTROPIC 
FLOW  EXPRESSIONS  WITH  FIXED  PROPESilES.  SMALL  SONIC  NOZZLES  HAVE 
MEASURED  DISCHARGE  COEFFICIENTS,  THE  AIR  FLOW  NOZZLE  USES  AN  ASSlRtED 
DISCHARGE  COEFFICIENT  (Cd)  OF  0.97, 

M  (LBM/SEC)=Cd»P*A»K«/T‘,5 

Kn  IS  THE  GAS-DEPENDENT  SONIC  CHOKE  FLOW  RATE  CONSTANT 
K«=SQR(  (Ganna*Gc/R)*(2/(Gahnatl))*((Gan'ia+l)/(Ganna-l))) 

APPROPRIATE  CONSTANTS  ARE: 


1030 

!  GAS 

MOLECULAR  WT 

GAS  CONST . 

CP 

GAMMA 

Kn 

1040 

1 

1050 

!  AIR 

28,97 

53.3 

.240 

1.40 

.5320 

1060 

!  02 

32,0 

48.3 

.217 

1.40 

.5589 

1070 

!  CH4 

16,03 

96.4 

.593 

1.32 

.3876 

1080 

!  N2 

28,01 

55,16 

.248 

1.40 

,5229 

1090 

!  H2 

2.016 

766.5 

3.419 

1.405 

.1405 

1100 

! 

1110 

Gannaair 

=1.40 

1120 

Ganna  01= 

1,40 

1130 

Gannahf= 

1.405 

Knair=.53EO 

Knho=.5589 
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1160  K«hf=.H05 
1170  Rair=53,3 
1180  Rho=4B.3 
1190  Rhf=766.5 
12D0  Dairrhiike=,42 
1210  Dbpchoke=.239 
12  Ohoctioke- .0700 
12j0  Dhfchoke=.040 
1241  Haird=l,9 
1250  Hbpaird=.59 
1260  Hfueld=.33 
1270  Hhfd=, 00285 
1280  Hhod=.0228 
1290  PRINT  USING  *6/* 

1300  INPUT  "Input  the  barohetric  pressure  in  nn  uf  Ho", Pun 

1310  Pbar=Pnn», 019337 

1320  Cdair=.97 

1330  Cdhf',97 

1340  Cdho=,97 

1350  !ALL  THERNOLOUPLES  ARE  CHRQHEL  us.  ALUMEL  (TYPE  K)  WITH 
1360  lELECTRONIC  ICE  POINTS.  TEMPERATURE  READINGS  (VOLTAGES)  ARE 
1370  ICONUERTED  TO  DEGREES  RANKINE  (R)  PER  "INDUSTRIAL  INSTRUMENTATION"  BY 
1380  !D.P.  ECKMAN  (PAGE  369),  THIS  CALCULATION  IS  PERFORMED  IN  SUBROUTINE 
1390  ITcalc,  TEN  VOLTAGE  INTERVALS  ARE  USED  BETWEEN  460  AND  2460  R, 

1400  PRINT  USING  "8" 

1410  PRINT  USING  "6/" 

1420  INPUT  "WILL  THE  AUGMENTQR  TUBE  BE  USED?  (Y/N>*,AuqO  '.DU 

1430  IF  Aufl$="Y"  THEN  Auq=l  iDU 

1440  IF  A;iq$=’N*  THEN  Aug=0 
1450  PRINT  USING  "8" 

1460  PRINT  USING  "6/* 

1470  INPUT  "WILL  THE  AIR  HEATER  BE  USED?  <Y/N>",Z2$ 

1480  IF  Zi$=’Y’  THEN  Ht=l. 

1490  IF  Z2$="N"  THEN  Ht=0. 

1500  PRINT  USING  "8" 

1510  PRINT  USING  *6/" 

1520  INPUT  "WILL  YOU  USE  PRE-INITIALIZED  VALUES  OF  CALIBRATION  CONSTANTS  AND  ZE 
R0S7  (Y/N)',Z2$ 

153d  IF  Zi<="Y"  THEN  GOTO  Initial 
1540  GOTO  Transcal 
1550  TcalC!  ! 

1570  i  »*»VOLTAGE  TO  TEMPERATURE  (RANKINE)  CONVERSION  SliPROUTINE»»i 

1590 
1600 
0+560 
1610 
7)+660 
1620 

2)  +760 
1630 

5) +860 
1640 
9)+960 
1650 

3) +1060 
1660 

4) +1160 
1670 

6) +1260 


IF  VoltsY, 00153  THEN  T=((Volts+. 00068)/. 00002201+460 
IF  Volt5)=. 00153  AND  Voltsl. 00382  THEN  T=( (Volts-. 00153)/. 000023 

IF  V<ilt5>=. 00382  AND  VoItsY. 00609  THEN  T=((Volts-, 00382)/. 800022 

IF  V(ilts>=, 00609  AND  Voltsl. 00831  THEN  T=((Volt5-. 00609)/. 000022 

IF  Vjlts>=, 00831  AND  Volts'.. 01056  THEN  T=((Volts-, 00831)/. 000022 

IF  Volts>=. 01056  AND  Volts(. 01285  THEN  T=((Volts-. 01056)/, 000022 

IF  Vi)lt5>=, 01285  AND  Voltsl. 01518  THEN  T=((Volts-. 01285)/. 000023 

IF  Volt5)=. 01518  AND  Vultsl. 01752  THEN  T=((Volts-. 01518)/, 000023 

IF  Volts>=, 01752  AND  Volts(. 01988  THEN  T=((Volts-. 01518)/, 000023 
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1680  IF  Voll5)=, 01988  AND  Volts<. 02225  THEN  T=l(yolts-. 01988)/, 000023 

7)+136D 

1690  If  Volts)=. 02225  AND  Vcltsl. 82463  THEN  T=((Vclts-. 02225)/, 000023 

B)+1460 

1700  IF  Volt5>=. 02463  AND  Volts(, 02698  THEN  T=((yoll5-, 02463)/, 000023 

5) +1560 

1710  IF  Volts)=, 02698  AND  Voltsl. 02932  THEN  T=((VoIfs-, 02698)/, 000023 

7) +1660 

1720  IF  y*lt8>=, 02932  AND  y»lts(, 03165  THEN  T=((Vb1is-, 02932)/, 000023 

3) +1760 

1730  IF  Volt5)=, 03165  AND  Voltsl. 03393  THEN  T=((yDlt5-, 03165)/, 000022 

8) +1860 

1740  IF  Vol1s)=. 03393  AND  VollsI. 03619  THEN  T=((Voll5-. 03393)/, 000022 

6) +1960 

1750  IF  Volt5>=, 03619  AND  Volls(, 03843  THEN  T=((Volts-, 03619)/, 000022 

4) +2060 

1760  IF  y#lls)=. 03843  AND  yoltsl. 04062  THEN  T=((yDlts-, 03843)/, 000021 

9)+2160 

1770  IF  yol1s)=, 04062  AND  Volts(, 04278  THEN  T=((Volts-, 04062)/, 000021 

6)‘*2260 

178o'’  IF  yiill5>=, 04278  AND  Voltsl. 04491  THEN  T=((Voll5-, 04278)/, 000021 

6) +2360 

1790  IF  yolts>=, 04491  THEN  T^dVolts-. 04278)/. 0000216)+2460 

1800  RETURN 
1810  Initial:  I 

1820  !  Initialued  values  of  leros  and  calibration  constants  for  all  transducers 
1830  Testno$='BEHRENS“ 

1840  D3te$=M'3fl-9fl* 

1850  Fuelidl=*0007" 

1860  Hedterfuel$='HYDROGEN' 

1870  Upa0=, 0017222 
1880  Kpa=33262,88 
1890  Vpc0=-, 176573 
1900  Kpc=844.814 
1910  ',Y)hi0=. 0356869 
1920  Kpho=67l. 61413 
1930  Vphf0=-.I16253 
1940  Kphf=1356,9902 
1950  Kfifuel  =  .5 
1960  Transcal;  ! 

1980  f(2)  TRANSDUCER  CALIPRAT10N3 

1990  'xiixxxxxxxxxxxxxxxxxxxxxoxxxxxxxixxtxxxtttxxxxxxxtxxxxxixafaxxxxxtxtxx 

2000  ! THERE  ARE  4  PRESSURE  TRANSDUCERS  THAT  hUST  BE  CALIBRATED 
2010  'TRANSDUCER  LINEARITY  MUST  BE  UERIFIED  BEFORE  THIS 
2020  ICALIBRATION  PROCEDURE  IS  EMPLOYED.  Tl€  ORDER  OF  CALIBRATION  IS  AS 
2030  IFOLLOyS;  Pa,  Pc,  Phf,  Pho 

2040  !THE  FOLLOUING  TWO  LINES  SET  UP  722  AND  709  FOR  DATA  ACQUISITION’ 

2050  ICl.EAR  709 
2060  CLEAR  722 
2070  REMOTE  709 

2080  OUTPUT  722;'LIRllSTNZ110STIT4gxI* 

2090  INPUT  ’DO  YOU  WANT  TO  CALIBRATE  TRANSDUCERS)  (Y/H)",Yy$ 

2100  IF  Yy$="N*  THEN  GOTO  Endcai 

2110  INPUT  ’DO  YOU  WANT  CONSECUTIVE  ORDER  (T  CALIBRAT10N7IY/N)'’,Yvt 

2120  IF  Yyt='Y"  THEN  GOTO  Consec 

2130  INPUT  ’DO  YOU  WANT  TO  RECALIBRATE  Pa)  (Y/N)",Yy$ 

2140  IF  Yy4="Y’  THEN  GOTO  Pacal 
2150  PRINT  USING 
2160  Pc:  I 
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2179 

2180 
2190 
2290 
2210 

a 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2339 

2340 
2350 
2360 
2370 
2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2560 
2570 
2580 
2590 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 


INPUT  'DO  YOU  WANT  TO  RECALIBRATE  Pc7  <Y/N)*,yy$ 

a»=’Y'  THEN  GOTO  Pccal 
USING  '8' 

IF  H1=0  THEN  GOTO  Endcal 
Phf;  ! 

INPUT  "DO  YOU  WANT  TO  RECALIBRATE  PhH  (Y/N>",Yyl 

IF  Yy9="Y'  THEN  GOTO  Phfcal 
PRINT  USING  "8" 

Phoi  ! 

INPUT  "DO  YOU  H.NT  TO  RECALIBRATE  Pho?  {y/N>",Yy» 

IF  Yy$="Y"  THEN  GOTO  Phocal 
GOTO  Endcal 
Consec:  ! 

Consul 
Pacal;  ! 

PRINT  USING  *2/' 

PRINT  "**  CALIBRATION  OF  Pa,  THE  AIR  SONIC  CHOKE  PRESSURE  TRANSDUCER*#" 
PRINT  USING  "2/" 

PaOcal:  ! 

PRINT  •  *****  ZERO  PRESSURE  «»*»*' 

PRINT  "INSURE  THAT  NO  PRESSURE  IS  APPLIED  TO  THE  TRANSDUCER" 

DISP  "HIT  CONTINUE  WHEN  READY  TO  TAKE  ZERO  READING" 

PAUSE 

REMOTE  799 
OUTPUT  709;"AC24" 

WAIT  2 

OUTPUT  722:  "Tl" 

ENTER  722;l)pa0 
PRINT  ■Upa0=‘;Upa0 
BEEP 

INPUT  'READING  OK?  (Y/Nl'.Zil 
IF  Z2t='N"  THEN  GOTO  PaOcal 
Paflaxcal;  ! 

PRINT  USING  "g" 

PRINT  '  »»»**  CALIBRATION  »*<*»" 

PRINT  'APPLY  MAXIMUM  PRESSURE  USING  THE  DEAD-HEIGHT  TESTER" 

INPUT  "ENTER  THE  MAXIMUM  PRESSURE  IN  psiq'.Pana* 

DISP  "HIT  CONTINUE  WHEN  READY' 

PAUSE 

REMOTE  709 
OUTPUT  709;'AC24' 

WAIT  2 

OUTPUT  722: "TS" 

ENTER  722)Vpa«a)( 

Print  "Vpana«=';Vp3Hax,'PaMa*=  "jPanax 
Kpa=(PaNax)/(Vp3nax-Vpa0) 

PRINT  ■Kpa=  'iKpa 
BEEP 

INPUT  'READING  OK?  <y/N)",Zzl 
IF  Zz»='N'  THEN  GOTO  Pafiaxcal 
IF  Consul  THEN  GOTO  Pccal 
Pccal;  ! 

I*********************************************************************** 

PRINT  USING  "8" 

PRINT  '**CALIBRATION  OF  Pc,  THE  T63  CHAMBER  PRESSURE  TRANSDUCER**" 

f*********************************************************************** 

Pc9cai;  ! 

PRINT  *«»»«  ZERO  PRESSURE  »*»»" 

PRINT  "INSURE  THAT  NO  PRESSURE  IS  APPLIED  TO  THE  TRANSDUCER" 

DISP  "HIT  CONTINUE  WHEN  READY" 
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I  PAUSE 
REMOTE  709 
I  OUTPUT  709 ; 'ACES' 

WAIT 

I  OUTPUT  722; 'T3' 

ENTER  722;')pc0 
I  PRINT  "VpcO='iVpcO 
BEEP 

INPUT  "READING  OR?  <Y/N>',Zl$ 

I  IF  Z2$='N'  THEN  GOTO  PcOcai 
Pcnaxcali  ! 

I  PRINT  USING 

PRINT  "  »«»»  CALIBRATION  *»*»' 

I  PRINT  "APPLY  THE  MAXIMUM  PRESSURE  USING  DEAD-WEIGHT  TESTER* 

INPUT  'ENTER  THE  HAXIMUH  PRESSURE  IN  psig'.Pcnax 
DISP  "HIT  CONTINUE  WHEN  READY' 

PAUSE 

REMOTE  709 
OUTPUT  709i'AC23' 

WAIT  2 

OUTPUT  722: 'T3" 

ENTER  722;VpcMax 

PRINT  •Vpc»iax=';Vpcflax,'Pc«ai=*;PcnaK 
Kcc=Pchax/(Upc«ax-Vpc(i) 

PRINT  ■Rpc=  ^Kpc 
BEEP 

INPUT  “READING  OK?  <Y/N)",Zii 
IF  Z2i=*N'  THEN  GOTO  Pcnaiical 
IF  Ht=0,  THEN  GOTO  Fincal 
IF  Cons=l  THEN  GOTO  Phfcal 
GOTO  Phf 
Phfcali ! 

PRINT  USING 

PRINT  'HCALIBRATION  OF  Phf  ,  THE  T63  AIR  HEATER  FUEL  TRANSDUCER  *»" 
Phfkal:  ! 

PRINT  'xxxxZCRO  PRESSURExxxx' 

PRiNT  'INSURE  THAT  NO  PRESSURE  IS  APPLIED  TO  THE  TRANSDUCER' 

DISP  'HIT  CONTINUE  WHEN  READY- 
PAUSE 

REMOTE  70? 

OUTPUT  709:"AC22" 

WAIT  2 

OUTPUT  722; 'T3' 

ENTER  722;')phf0 
PRINT  'Vphffl=';VphfO 
BEEP 

INPUT  "READING  OK?  (Y/H)',Zz$ 

IF  Zz»='N‘  THEN  GOTO  Phf Oca  1 
Phfflaxcal:  ! 

PRINT  USING  ’P" 

PRINT  "xxxxCALIBRATIONxxxx* 

DISP  'APPLY  THE  MAXIMUM  PRESSURE  USING  DEAD-WEIGHT  TESTER' 

INPUT  'ENTER  THE  MAXIMUM  PRESSURE  IN  psio'.Phfnax 
DISP  'HIT  CONTINUE  WHEN  READY' 

PAUSE 


3360  REMOTE  709 
3370  OUTPUT  70?)'AC22' 
3380  WAIT  2 
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3390  OUTPUT  722: 'TS* 

3408  ENTER  722iVphf«ax 

3411  PRINT  ■Vphfha*=":Vphf«ax  "PhfM«='}PhfHax 
3420  KphfrPhfnai/lUphfMX-Vphh) 

3430  PRINT  •Kphf=*il(phf 
3440  PEEP 

3450  INPUT  'READING  OR?  (Y/N)",Z2» 

3460  IF  22$=‘N’  then  GOTO  PhTnaical 
3470  IF  Con5=l  THEN  GOTO  Phocal 
3480  GOTO  Pho 
3490  Phocal;! 

3510  PRINT  USING  *0" 

3520  PRINT  'ixCALIBRATION  OF  Pho,  THE  AIR  HEATER  OXYGEN  PRESSURE  TRANSDUCER »*• 

3530  IXXXXXXXMXXXIXKXVXXXXtXIxilMIMXXXXXtRXHKXfXXXXXKXHXXXIXXHtXXKIXIXX’X 

3540  Ptioflcal:  ! 

3550  PRINT  'xxxxZERO  CALIPRATION****’ 

3560  PRINT  'INSURE  THAT  NO  PRESSURE  IS  APPLID  TO  THE  TRANSDUCER' 

3570  DISP  "HIT  CONTINUE  UHEN  READY' 

3580  PAUSE 

3590  REMOTE  709 

3600  OUTPUT  709; 'ACPI' 

3618  OUTPUT  722; "T3' 

3620  ENTER  722)Vpho0 

3630  PRINT  'Vpholl=*;VphoO 

3540  INPUT  "READING  OK?  (Y/N)-,Zz$ 

3650  IF  ZzS^'N'  THEN  GOTO  PhoOcal 
3660  Phonaxcal;  I 
3670  PRINT  USING  '§' 

3680  PRINT  “»»«»CALIBRATIOH*«»*' 

3690  PRINT  'APPLY  THE  HAXINUM  PRESSURE  USING  DEAD-WEIGHT  TESTER' 

3700  INPUT  ’ENTER  THE  NAXIHUM  PRESSURE  IN  psig’,PhoNax 

3710  DISP  'HIT  CONTINUE  WHEN  READY' 

3720  PAUSE 
3730  REMOTE  709 
3748  OUlPUT  709;*AC2r 
3750  OUTPUT  722; "T3' 

3768  ENTER  722;0pho«ax 

3770  PRINT  'Uphonax=";'7phonax,'Phonax='jPhoMax 
3780  Kpho=Pho«ax/(UphoMX-VphoO) 

3790  print  'Kpho^'iKpho 
3800  PEEP 

3810  INPUT  "READING  OK?  {Y/N)",Z2$ 

3820  IF  Z2»="N’  THEN  GOTO  Phonaxcal 
3830  IF  Auq=l  THEN  GOTO  Auqcal 
3840EndcaI;  ! 

3850  Fincal;  I 
3860  PRINT  USING  "8" 

3870  PRINT  'THIS  ENDS  THE  CALIBRATIONS  " 

3890  !(J)  PRE-RUN  INPUTS,  FLOW  RATE  CHECKS  AND  NOZZLE  CALCULATIONS 

3909  !XXX»XXX»>X«X>XXXXXMXXXXi»»X»»>XX»XXXIXXX«»«XiXX>IXXX»XX(«IXX«IXXXX'« 

3910  IXXIKXtXXXIXXtXXXfXtXXXXXXXXXXXXXXXXXXiXHXXXXXXXXXXffXXXXXXXiXXXXXXfiXXXX 


3920  !  A,  FLOW  RATE  SET-UPS  AND  CHECKS 
3940  PRINT  USING 

3950  INPUT  'DO  YOU  WANT  TO  PRESET  THE  AIR  FLOW  RATE?a/N>',Z2$ 

3960  IF  Z2»="N'  then  GOTO  Paskip 

3970  PRINT  'SET  THE  DESIRED  VALUE  OF  Palpsig)  USING  THE  HAND  LOADER  /PRESSURE  G 
AGE' 

3980  PRINT  'USING  '3/' 

3990  PRINT  'THE  HAND  LOADER  SHOULD  PE  20  PSIG  MORE  THAN  DESIRED  PRESSURE' 
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Paset!  ! 

PRINT  USING  'ir 

PRINT  "MANUALLY  INITIATE  AIR  FLON  BY  TURNING  'MAIN  AIR'  ON  CONTROL  PANEL' 
PRINT  USING  *2/" 

DISP  "HIT  CONTINUE  WHEN  READY" 

PAUSE 
WAIT  3 

OUTPUT  709j*AC24" 

OUTPUT  722; "TS" 

ENTER  722; Vp a 
OUTPUT  70^/ACbO" 

OUTPUT  722; ’13 " 

ENTER  722iUta 
OUTPUT  709; "ACS r 
OUTPUT  722; 'T3" 

ENTER  722;(lti 
BEEP 

PRINT  USING  "e* 

PRINT  "TURN  OFF  'MAIN  AIR'" 

DISP  "HIT  CONTINUE  TO  FRIKEED" 

PAUSE 

Pa=(Upa-VpaO)»Kpa+Pbar 

V«lts=Uta 

GOSUB  Teak 

Ta=T 

Volts=Vti 
GOSUB  Teak 
Tcin=T 

Mair=Knair»Cdair»Pa*.785A*(Dairchoke*2)/(Ta* .5) 
Ppair=Knair»Cdair*Pa».7854*(Dbpchoke*2)/(Ta''.5) 

PRINT  USING  "8" 

PRINT  USING  "5A,2X.DDD,DDDD';"Mair=*;Hair 
PRINT  USING  •14A,DDD,DDD";"Hair  DESIfiED=";«aird 
Ratii)=Mair/Maird 

PRINT  USING  •20A,D,DDD,2X,3A,lX,DDDD.D,lA,3X,3A";'Mair/DESIRED  HaiP=";Rat 
Ta='‘Ta  ‘R*  »  i  <  < 

P((lNf  USING  •6A,2X,DDD.DDDD';'BPair=*)Bpair 
PgsPa-Pbar 

PRINT  USING  "4A,DD0D,D,6A";'Pa  =*;Pg;"  Psiq* 

PRINT  USING  "4A,DDDD,D.3A"j"Ta  V 

INPUT  "IS  AIR  FLOW  RATE  ACCURATE  ENOUGH?  (Y/N)',Xx$ 

IF  Xx$="Y'  THEN  GOTO  Prerun 
Paneu=(Pai(Maird/Mair)-Pbar 
PRINT  "RESET  Pa  TO’;FaneMj'PsiQ" 

DISP  'HIT  CONTINUE  AFTER  RESET 'OF  Pa" 

PAUSE 

GOTO  Paset 
Prerun:  ! 

INPUT  'DO  YOU  WANT  PRINTOUT  OF  PRE-RUN  DATA?<Y/N)",Xx» 

IF  Xx$="Y"  THEN  GOTO  Preprint 
GOTO  Skipprint 
Preprint;  ! 

PRINTER  IS  701 

PRINT  •  *»»*  PRE-RUN  DATA,  USING  AIR  ONLY  »***" 

PRINT  " 

PRINT  'DATE:  ";Oale$ 

PRINT  USING  •3A,DDDD,D,6A"j’Pa="iPa,"  Psia’ 

IPRINT  USING  '3A,DD0,0,3A'|'Tcin=%Tcin,’  R" 

PRINT  USING  '5A,D,DD0D,llll";"Hair=’:Hair,’  ab«/sec)* 

PRINT  USING  ’5A,D.DDDD,llA'i"BPair=';Bpair,’  (Ibn/sec)" 
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4590  PRINTER  IS  I  . . 

4680  Skipprint;  ! 

4610  DISP  'HIT  CONTINUE  TO  PROCEED  TO  NEXT  FLON  RATE  SET  UP* 

4620  PAUSE 
4630  Paskip:  ! 

4640  IF  Ht=0  THEN  GOTO  Phoskip 

4650  PRINT  USING  ’S* 

4660  INPUT  'DO  YOU  WANT  TO  PRESET  THE  HEATER  FUEL  FLOH  RATE?  {Y/N>*,Z2» 

4670  IF  Zzl=*N*  THEN  GOTO  Phfskip 

4690  PRINT  'SET  THE  DESIRED  VALUE  OF  Phf  USING  THE  HAND  LOADER/PRESSURE  GAGE" 

4710  DISP  ‘HIT  CONTINUE  HHEN  READY" 

4720  PAUSE 
4730  Phfset : ! 

4740  PRINT  USING  -?* 

4750  PRINT  'MANUALLY  TURN  ON  AIR  'HEATER  FUEL'  SNITCH" 

4760  DISP  ■  HIT  CONTINUE  TO  PROCEED* 

4770  PAUSE 

4780  OUTPUT  709;*AC22" 

4790  OUTPUT  722; *T3" 

4800  ENTER  722:0phf 
4810  OUTPUT  709;'aC65* 

4620  OUTPUT  722: 'TI" 

4830  ENTER  722;Vthf 

4840  ICLEAR  709 

4850  PRINT  ‘MANUALLY  TURN  OFF  AIR  'HEATER  FUEL'  SNITCH* 

4860  PEtP 

4870  DISP  'HIT  CONTINUE  TO  PROCEED* 

4880  PAUSE 

4890  Phf=<Vphf-VphfO)»Xphf+Pbar 
4900  Volts=Vthf 
4910  GOSUB  Teak 
4920  Thf=T 

4930  Hhf=K«hf»Cdhf*Phf*.7854*(Dhfchoke'‘2)/(Thf‘.5) 

4940  PRINT  USING  ‘8* 

4950  PRINT  USING  "4A,DD,DDDDD";*Hhf=";Mhf 

4960  PRINT  USING  'IBA.DD.DDDDDV 'Mhf  flESIRED=';«hfd 

4970  Raijo=Mhf/Mhfd 

4980  PRINT  USING  *18A,D.DDD,2X,4A,3DDD,DD,lA’|"Mhf/  Mhf  DESIRED="iRatiii,*Thf=* 
499o’  *Pg=Phf-Pbar 

5000  PRINT  USING  "5A,DDDD.DDD,4A,3X,4A,DDDD.DD,lA,4A,DDDD,DD,lA*;*Fhf=  ’iPqi'P 
sig',"Thf="jThfj’R" 

50ffl  INPUT  "IS  HEATER  FUEL  FLOH  RATE  ACCURATE  ENOUGH?  a/N>*,Xx» 

5020  IF  Xx4=*Y*  THEN  GOTO  PhPfin 
5030  PhfneM=(Phf*Mhfd/Hhf>-?bar 

5040  PRINT  USING  "I3A,DPPD.DD,4A": ‘RESET  Phf  TO";Phfnew;'P5ig’ 

5050  DISP  ‘HIT  CONTINUE  AFTER  RESET  OF  Phf 
5060  PAUSE 
5070  GOTO  Phfset 
5080  Phffin;! 

5090  DISP  'HIT  CONTINUE  TO  PROCEED  TO  NEXT  FLOH  RATE  SET  UP" 

5100  PAUSE 

5110  PhfskiD;  ! 

5120  PRINT  USING  ‘8* 

5130  INPUT  "DO  YOU  HANT  TO  PRESET  THE  HEATER  OXYGEN  FLOH  RATE?(Y/N>' ,Zi$ 

5140  IF  Zz»="N’  THEN  GOTO  Fboskip 

5160  PRINT  ‘SET  THE  DESIRED  VALUE  OF  Pho  USING  THE  HAND  LOADER/PRESSURE  GAGE" 

5170  IXXXiXXKXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXiXXXXXXXXiXKXXXXXXXKXXXXXtXXXXXX 

5180  Phoset:! 

5190  PRINT  ‘MANUALLY  TURN  ON  AIR  'HEATER'  OXYGEN'  SHITCH* 

5200  DISP  ‘HIT  CONTINUE  TO  PROCEED* 
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PAUSE 

OUTPUT  709; 'AKl' 

OUTPUT  722; "TS* 

ENTER  722;(ipho 
OUTPUT  709;'aC6A* 

OUTPUT  722!’T3'' 

ENTER  722;Utho 

PRINT  'HANUALLY  TURN  OFF  AIR  'HEATER  OXYGEN'  SWITCH" 

BEEP 

DISP  "HIT  CONTINUE  TO  PROCEED* 

PAUSE 

Pho=(Vpho-VphoO)*Kpho+Pbar 

Vfllts=Vtho 

G05UB  Teak 

Tho=T 

Hho=K«ho*Cdho*Pho«.7854»(Dhochoke*2)/(Tho“.5) 

PRINT  USING  "8" 

PRINT  USING  '4A,DD,DDDDD":*Hho=":ttho 
PRINT  USING  •18A,DD.DDDDD'';”Mho  DESIRED=*;Hhod 
Batio=Nhfl/Hhod 
Pq=Pho-Pbar 

PRINT  USING  *5A,DDDD.DD,tX,5A,5X,4A,DDDD.DD,tX,2A*i"Pho=*iPq;"Psig’i 


"Tho=';Tho;‘R"  >  ■  • 

5430  INPUT  "IS  THE  HEATER  OXYGEN  FLOW  RATE  ENOUGH?  (Y/N>?",Xx$ 

5440  IF  X)(»="Y"  THEN  GOTO  Phoskip 

5450  PhoneM=(Pha*Mhod/Mho)-Pbar 

5460  PRINT  USING  "14A,DDDD,DD,}X.4A*: "RESET  PhD  TO  ")PhDnen;'Psig’ 

5470  DISP  "HIT  CONTINUE  AFTER  RESET  OF  Pho" 

5480  PAUSE 

5490  GOTO  Phoset 

5500  Phoskip !  ! 

5510  PRINT  "THIS  COMPLETES  PRE-RUN  SET-UP* 

5530  !  <4)  THIS  PORTION  OF  THE  PROGRAM  RUNS  THE  TEST  AND  COLLECTS  THE  DATA 

5540  jXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXIXXXIXXXXXXXXXX'XX 


5550  !  PRINT  USING  "8" 

5560  DISP  "SET  TIMEDATE  BY  PRESSING  K19  AND  UPDATE;  THEN  EXECUTE,  THEN  HIT  CON 
TINUE" 

5570  BEEF 

5580  PAUSE 

5590  PRINT  USING  "8* 

5600  Rpt!  ! 

5610  PRINTER  IS  I 

5620  !  THE  FOLLOWING  PROGRAMS  THE  3456  DVM 

5630  ASSIGN  SScanner  TO  709 

5640  ASSIGN  8Svh  TO  722 

5650  CLEAR  8Svn 

5660  OPTION  BASE  1 

5670  DIM  Press! 10,5) 

5680  DIM  TenpdO.t) 

5690  DIM  Auqfeflp(l0,3) 

5700  CLEAR 

5710  OUTPUT  8SvM;"LlZlD0S0FlR30STDlSTNP0rL01STISOlT401QXl" 

5720  DISP  "HIT  CONTINUE  FOR  HOT  RUN  DATA" 

5730  BEEP 

5740  PAUSE 

5750  DISP 

5760  GOSUP  Press 

5770  GOSUB  Ten 

5780  IF  Aug=l  THEN  GOSUP  Aug 

5790  GOTO  shutdown 
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5800 

5810  Press;  ! 

5830  PRINT  USING  "8" 

5830  PRINT  USING  MO/" 

5840  PRINT  ’  COLLECTING  PRESSURE 

5850  PRINT  ■■ 

5860  PRINT  •  »*»*»«»*»»**  COLLECTING  PRESSURE  »»««»»»*»»«»’ 

5670  OUTPUT  8Scanner;-AC31AF21AL25AE2“ 

5880  WAIT  .2 

5690  OUTPUT  8Svn!*50STNT3" 

5900  ENTER  eSv«  USING  "♦.K'jPressI*) 

5910  !  SUBROUTINE  Press 

5920  OUTPUT  8Sv«;"lSTNT4' 

5930  RETURN 

5940  !Mtt*4Mt0t444ttH4O44»0H4Hm4«OtHHH4«4«04H4**«HHtHm*M»44O 

5950  Aug:  !  TEMPERATURE  COLLECTING  ROUTINE  FOR  AUGHENTDR  TUBE 

5960  PRINT  "M*  COLLECTING  AUGMENTOR  TEMPERATURES  ««*»■ 

5990  OUTPUT  RScannerj •AC66AF66AL6BAE2’ 

5991  OUTPUT  8Sg«;"30STHR2T3" 

5992  «AIT  .2 

6000  ENTER  eSvn  USING  '».K";Auqtenp(») 

6010  OUTPUT  eSvft)  "RSI STNT4* 

6020  RETURN 

6030  !«*44444MfHM444444Mf44MH4H40<i«Mtt444*444MH«*t*»4»4t4*MtM4t 

6040  Ten;  ! 

6050  PRINT  •  *»»*»COLLECTING  TEMPERATURES  *««»*«- 

6060  OUTPUT  8S»n;'60STNR2" 

6070  OUTPUT  8Scanner;*AC60AF60AL65AE2" 

6080  NAIT  ,2 

6090  OUTPUT  8Svii;’60STNR2T3" 

6100  ENTER  8Sv«  USING  ’'l.i:" iTeMpi*) 

6110  OUTPUT  8Sv«;’lSTNR3T4‘ 

6120  RETURN 

6130  !44444«Mt4tH*M44444M44tMItM44ttHMM»»4t44H0HMM«tH4H444t44t 

6140  Shutdtun;  ! 

6150  PRINT  USING  "8" 

6160  PRINT  '  TEST  COMPLETE;  TURN  OFF  MAIN-AIR,  HEATER  GASES  ' 

6180  !»  (5)  POST-RUN  OPERATION,  DATA  REDUCTION  AND  SHUTDOWN 

6190  ixxxxxxxxxxxxxxxxxtxxxxxxyiixyyxxxyxxxxxxxxxKXxixX'HxxxixxtHXXHiX'XtxiHxxx 

6200  DISP  "HIT  CONTINUE  TO  PROCEED  TO  DATA  REDUCTION" 

6210  BEEP 

6220  PAUSE 

6230  PRINTER  IS  701 

6240  PRINT  USING  '3/' 

6250  PRINT  '■  «***  PRE-RUN  INPUT  »*»»  ’ 

6260  PRINT  USING  ’2/" 

6270  PRINT  USING  ■14A,?A,5X,14A,9A,5X,14A,9A';'Te5lno=";Te5tno$,’Dale="iDate»,' 
Fuelid=";Fuelidl 

6280  PRINT  USING  ’14A,9A" :''Heaterfuel=";He3terf«el» 

6290  PRINT  USING  "14A,DDDD,DDD,6X,14A,DDDD,DDD";"Dairchcl(e=";Dairchclie,"Dhtfuch 
oke=";Dhfchol(e 

6300  PRINT  USING  "14A,DDDD,DDDD':"Dhtoxch()l(e=*;Dhochoke 

6310  PRINT  USING  "14A,DDDD.DDD,6)(,14A,DDDD,DDD  7X,14A,DDDD,DDD"; •Cdair=";Cdair, 

"Cdh=*;Cdhf,’Cdho='iCdf!o 

6320  Print  using  •i4A,DDDD.DDD,5X,14A,DDDD.DDD,5X,14A,DDDD.DDD’;"GanHahf=";Gahn 
ahf 

6330  PRINT  USING  *14A,DDDD.DD.5X.14A,DDDD,DD.5X.14A,DDDD.DD’;'Rhf="iRhf 
6340  PRINT  USING  '14A,DDDD,DIDD,5X,14A,DDDD.DDDD,5X,14A,DDDD.DDDD";''l(Mhf="l*nhf 
6350  PRINT  USING  "HA.DDDD.DDDD  5X:i4A.DDDD,DDDDM'l(«ho=":Kflho.’KHair=":KHair 
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360  PRINT  USING  "}4ft,DDDD.DDD.iX.\4A.t)DDDp*:’Naird=;!Haird, "tcind=';Tcind 
370  PRINT  USING  M4A,DDDD.DDDi),5ll,14^.DDDD.DllDD,5X,14A,DDDD,DDDD,4X,I4A,DDDD.D 
iDD’ : ’Hhfds"  i«hf  d ,  ■Hhod=' iMhod . ’nf  ueld=^Hf  ueld 

dsiNG  ■  14A,D0iDD . d4i>D,5X , \ 4A  DDMD .OMD  .5X  A 4A ,DDDDD .CDDD ,5X ,  1 4A,DDD 
a=*^l(ga^*VgaO=*jVpaO;'ipc=*jKpcj*Vpc«=-;Vpcll 
^  “  3/ 


6 
6 

DDD 
6380 

D.DDDD'j*Kpa=*jKp 


63?0 

6406 

6411 

6420 

6430 


PRINT  USING 
PRINT  " 

PRINT  USING  ■35X,28A'j'«h»  DATA  EXTRACTED  «*»** 

PRINT  •" 

PRINT  USING  M2DA'i’FUEL  FLOU  RATE  IN  Cal/nin,  GAS  FLOU  RATES  IN  Lb«/sec, 
PRESSURE  IN  Psia,  TEMPERATURES  IN  R' 

6440  PRINT  USING  ■5A.6X.6A,5X,6A.7X.6A,9X,6(7A,5X)’;’Count',*Mair’,'«bpair 
■Hfoer/f/Hhf*,‘«ho*/Phr,‘Pho',‘Pc* 

'■  FdR  J=1  TO  fo 
Vpa=PressU.4) 

Pa=(Vpa-VpaO)»Kpa+Pbar 
Vpc=Pre5s(J,3) 

Pc=(Vpc-Vpcu)iKpc+Pbar 
Vnftfel=Pre5s(j,5) 
lifgel=Uiifuel»l(«fuel 
IF  H1=1  THEN 
Vphf=Press(J,2) 
piif=(Vphf-VphfO)»Kphf+Pbar 
Vpho=Pres5U,l) 

Ph«=(Vpho-VphaO)»Kph#+Pbar 
END  IF 

Vol1s=Te«p(J,l ) 

GOSUB  Tcalc 
Ta=T 

IF  HIM  THEN 
VoIf5=Te«p(J,5) 

GOSUB  Tcalc 
ThcM 

VoIts=Teflp(J,6) 

GOSUB  Teak 
ThfM 
END  IF 


6^50 

6468 

6470 

6480 

6470 

6500 

6510 

6520 

6530 

6540 

6550 

6560 

6570 

6580 

6590 

6600 

6610 

6620 

6630 

6640 

6650 

6660 

6670 

6680 

6690 

6700 

6710 

6720 

6730 

6740 

6750 

6760 

6770 


Mair=Kfiair»Cdair<Fa«.7854<(Daircholie‘2)/{Ta\5) 

Bpair=K«air*Cdair*Pa».7854*(Dbpchoke‘2)/(Ta*.5) 

F=,108»Mf«el/Hair 

IF  Ht=e,  THEN  GOTO  Junp 

IF  Phf(MPc»2)  THEN  Hhf=0. 

IF  Fho<MPc«2)  THEN  tIhoM , 

Hhf=Knhf«Cdhf«Phf»,7854»(Dhfch(ike*2)/(Thf‘.5) 

Hho=Knho*Cdho*Pho*.7854»(Dhochoke‘2)/(Tho'.5) 


6780  Junp;  ! 

6790  IF  H1=0,  THEN 

6000  Phf=0 

6810  Pho=0 

6820  Hhf=0 

6830  ltt)C=0 

6840  END  IF 

6850  !  THE  VALUES  BELON  ARE  DEFAULT  VALUES  TO  PREVENT  PRINTER  ERRORS 

6860  IF  Nairl.OOOl  THEN  Nair^.OOOl 

6870  IF  Bpairl.flOOl  THEN  Bpair=,0001 

6880  IF  MfueU.OBOl  THEN  Hfvek.OOOl 

6890  IF  F<. 0001  THEN  F=. 0001 

6900  IF  «hf(. 0001  THEN  Hhf=, 0001 

6710  IF  Mhol. 0001  THEN  Hho=. 0001 

6920  IF  Phf<, 0061  THEN  Phf=,0001 

6930  IF  Pht<. 0001  THEN  Pho=. 0001 

6940  IF  PcI.Oflfll  THEN  Pc=. 0601 
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6950  IF  Hair)2000,  THEN  t1air=2000. 

6960  IF  BDair)2000.  THEN  Bpair=2DD0. 

6970  IF  Mruel>2800.  THEN  Hfuel=2000. 


6980  IF  F>2000,  THE*  F=2000, 

6990  IF  «hf>2fl0fl.  THEN  Mhf=2fl(10, 

7000  IF  Hho)2000.  THEN  Hhf=2000. 

7010  IF  Phf >2000.  THEN  Phf=2000, 

7020  IF  Pho)200fl.  THEN  Pho=20(l0. 

7030  IF  Pc>200fl.  THEN  Pc=2000, 

7040  PRINT  USING  "DDD,3X,HD.DDDE  JX.HD.DBDE  ,3X,I1D.DDDE  ,3X,3(HD.DDDE  ,3X 
).3(«DD.DDE  ,3X),lX*jJ,Hair,Bpair,Hruei,F,Hhf,Mh#,Phf,Pho,Pc 
7(l50  NEXT  J 

706(  PRINT  "• 

7070  PRINT  USING  ■5A,5X,9(7A,7X)';"Counf,*TaV'Tcin',’Texr,'Tex2'  ’Taugup 

*.'Taugdr,'TaBqd2" 

7080  FOR  J=rTO  10 

7090  Volts=Te«p(J,l) 

7100  GQSUB  Tcalc 

7110  Ta=T 

7120  Uolts=Tehp(J,2) 

7130  GOSUB  Tcalc 

7140  Tcin=T 

7150  Volts=Te«p(J,3) 

7160  GQSUB  Tcalc 

7170  Te»l=T 

7180  Uolt5=Te«p(J,4) 

7190  GOSUB  Tcalc 

7200  Tex2=T 

7210  Volt5=Augte«p(J,l) 

7220  GOSUB  Tcalc 

7230  Taugup=T 

7240  Volts=Auqte«p(J,2) 

7258  GOSUB  Tcalc 

7260  Tauqdl=T 

7270  V(ilts=AuqteMp(J,3) 

7280  GOSUB  Tcalc 

7298  Tai>gd2=T 

7300  IF  Ta<100.  THEN  Ta=100. 

7310  IF  TcinllOO,  THEN  Tcin=100, 

7320  IF  TexlUOO,  THEN  Texl=100. 

7330  IF  Tex2(10C,  THEN  Tei2=100, 

7340  IF  Taug#p<l03.  THEN  TaugBp=100. 

7350  IF  TaugdHlOO,  THEN  Tauadl=180. 

7360  IF  Taugd2<10fl.  THEN  Taudd2=t0fl. 

7370  IF  Ta>408fl.  THEN  Ta=40C0, 

7380  IF  Tcin>4000,  THEN  Tcin=4000, 

7370  IF  Te:l)4C00.  THEN  TsxlMCOO. 

7400  IF  Tex2>4flOO,  THEN  Tex2=4000. 

7410  IF  TaBgup>4000.  THEN  Ta*gup=4000. 

7420  IF  Taugdl>4000.  THEN  TaugdlMOOO. 

7430  IF  Taugd2)4008.  THEN  Tai/gd2MPffl. 

7440  PRINT  USING  ■DDD,3X,7(MDD.DDE,4X)";J,Ta,Tcin,Texl,Tex2,Taugijp,Taugdl,T 

auqd2 

7450  NEXT  J 

7460  GOTO  Finish 

7470  Finish:  ! 

7488  PRINTER  IS  1 

7490  PRINT  USING  '8" 

75(0  PRINT  'DATA  OUTPUT  IS  COMPLETE’ 

7510  DISP  ’SECURE  TEST  CELL  !!!’ 

7520  BEEP 

7530  PAUSE 

7540  END 
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APPENDIX  B 


RUN  CHECKLIST 


TEST  CELL  #1 

1.  Ensure  yellow  and  top  blue  air  valves  in  the  solid  fuel 
ramjet  test  cell  are  closed. 

2.  Open  lower  blue  valve  (opens  air  line  to  Test  Cell  #2  or 
T-63)  . 

**  Note  -  At  least  one  valve  should  be  open  at  all  times  from 
the  main  air  line  to  ensure  an  air  vent  in  case  of  an 
accidental  component  failure. 


NITROGEN  BOTTLE  ROOM 

1.  Fully  open  the  control  room  nitrogen  bottle.  Ensure  that 
there  is  at  least  1000  psi  available. 

2.  Fully  open  actuator  nitrogen  bottle.  Ensure  that  there 
is  at  least  500  psi  available. 


CONTROL  ROOM 

1.  Ensure  AC  master  switch  is  on  and  the  red  covered  main  air 
switch  is  closed  on  the  solid  fuel  ramjet  control  panel. 

2 .  Ensure  there  is  no  pressure  set  on  the  gauge  on  the  Air 
Heater  console  (the  air  flow  set  pressure) . 

3 .  Ensure  that  the  T-63  combustion  chamber  safety 
thermocouple  is  installed  and  operating. 

4.  Ensure  that  the  fuel  tank  set  pressure  (gauge  on  panel) 
is  less  than  500  psi. 


FUEL  STORAGE  ROOM 

1.  Open  nitrogen  bottle  valve  (need  at  least  400  psi  more 
pressure  available  in  the  bottle  than  the  desired  fuel 
line/tank  pressure,  or  900  psi  minimum) . 
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2.  Adjust  hand  loader  to  read  700  psi. 

3.  Slowly  open  the  nitrogen  gas  supply  valve  located  behind 
the  fuel  tank  near  the  wall. 

4.  Very  slowly  open  the  fuel  line  valve  from  near  the  bottom 
of  the  tank  to  the  T-63. 


OUTSIDE/CONTROL  ROOM 

1.  Open  main  air  valve  to  full  open  (ensure  that  there  is  at 
least  2500  psi  available  in  outdoor  air  storage  tanks  for  a 
run. 


2.  Ensure  all  thermocouples  are  turned  on  (if  required)  and 
pressure  transducers  and  tubing  are  secure  at  the  test  stand. 

3.  The  heated  sample  line  temperature  control  box  should  be 
set  to  275  deg.  F  and  the  gas  analyzers  in  the  control  room 
should  be  up  and  operating.  The  three  main  switches  for  the 
electronic  equipment  racks  should  be  ON. 

4.  Load  and  run  the  "T63NOX"  computer  program  on  the  HP 
microcomputer.  The  pressure  transducers  should  be  calibrated 
if  not  already  done  so,  and  zeros  and  constants  entered. 

5.  Set  the  main  air  pressure  to  600  psi  using  the  hand 
loader. 

6.  Set  the  fuel  tank  pressure  to  500  psi  using  the  hand 
loader. 

7 .  Go  through  the  flow  rate  set  procedures  in  accordance  with 
the  computer  program.  This  may  require  opening  and  resetting 
the  pressures  for  the  air  heater  fuel  and  torch  gases  as  well 
as  the  heater  oxygen  tanks  in  the  test  cell. 

8.  Ensure  printer  is  "on-line”. 

9.  Check  for  personnel  near  the  test  cell  and  for  golfers. 
Activate  exterior  warning  horn  and  check  main  air  flow  rate 
as  cued  by  the  computer. 

10.  Turn  on  rocket  motor  siren. 

11.  Start  strip  chart  recorder  and  mark  zero/ambient 
conditions. 

12.  Activate  main  air  ON. 
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13 .  Turn  on  air  heater  fuel  and  heater  torch  (momentary)  until 
light-off,  if  required. 

14.  Ensure  T-63  engine  ignitor  key  is  installed  and  in  the  ON 
position. 

15.  Simultaneously,  activate  the  toggled  engine  ignitor  switch 
and  fuel  switch.  Check  desired  fuel  flow  rate  (0.33  GPM)  . 
Watch  for  wet  or  hot  start  by  visually  observing  exhaust  smoke 
at  rig  and  monitoring  the  digital  combustion  chamber  safety 
temperature  readout  (commence  shutdown  if  temperature  reaches 
1380  deg.  F) 

16  When  steady-state  operation  is  reached,  begin  traversing 
the  Kiel  probe  in  the  augmentor  tube  and  obtain  analyzer 
measurements . 

17.  After  data  gathered,  switch  fuel  OFF  and  air  heater  OFF 
(if  applicable) .  Leave  main  air  ON  until  engine  and  augmentor 
tube  are  cool. 

18.  Turn  main  air  OFF,  record  run  time,  and  calculate  fiiel 
used  during  run.  Update  fuel  board  in  fuel  storage  room. 

19.  Isolate  fuel  tank  with  valves  and  bleed  excess  fuel  in 
lines  with  fuel  switch  activation. 

20.  Secure  all  torch  and  air  heater  gas  bottles  in  test  cell. 

21.  Close  main  air  valve  outside. 

22.  Vent  fuel  tank  from  control  panel  if  desired  and  close 
fuel  tank  nitrogen  bottle. 

23.  Bleed  remaining  air  heater  and  torch  gases  from  lines  and 
vent  with  remaining  main  air  in  lines.  Back  off  pressure 
loaders  to  zero  in  the  control  room. 

24.  Secure  analyzers,  complete  shutdown,  and  reduce  data. 
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APPENDIX  C 


HOT  RUN  DA-^A 

DATE;  3-|f-9il 
Pa=  573.1  Psia 
Mair=1.9048  (Lbti/sec) 

BPair  .6168  (Ibn/sec)  _  i^q  heater 

No  catalyst  used 


«***  PRE-RUN  INPUT  «*«* 


Testno- 

BEHRENS 

Date= 

S-l4f-9fl 

Fuelid^ 

Heaterfuel= 

HYDROGEN 

Dairchoke= 

.420 

DhtPuchokg: 

.040 

Dhtoxcholie= 

.0700 

Cdair= 

Ga(4(tahf= 

.970 

1,405 

Cdh= 

.970 

Cdho= 

Rhf= 

766,50 

Ktihf= 

.1405 

Knho= 

.5589 

K«air= 

,5320 

Maird= 

1.900 

Tcind= 

0,0 

llhfd= 

,0028 

Mhod= 

.0228 

Hfueld= 

Kpa= 

33262,8800 

Vpai,= 

.0017 

Kpc= 

****  DATA  EXTRACTED  «»** 


1007 


.970 


.3300 

844.8140 


Pc 

92,<3E+00 


92.11Eti)0 

92.62E+00 


92,28E+00 

92.J8E+00 

92.35E+00 

91.24E^D0 


92.05E400 

91.62E40D 

91,42E+00 


FUEL  FLOW  RATE  IN 


Count 

1 

n 

L. 

3 

4 

5 

6 

7 

8 
9 

10 


Nair 
1. 9226400 
1.9l9Et20 
1.912E4C0 
1.9l3E+fl0 
1.923E400 
1.918E400 
1,920E400 
l,919E+00 
I,920£+00 
l,S22E+00 


Gal/Min,  GAS 
Mbpair 
6,224£-01 
6,216E-01 
6,I9I£-0I 
6.193E-01 
6,225E-0; 
6.211E-01 
6.216E-01 
6,214E-01 
6,2i6E-01 
6.222E-01 


FLOW  RATES  IN 
Nfuel 
3.026E-01 
3.C39E-01 
3,043E-01 
3.fl44E-fll 
3.fll5E-01 
3.045E-01 
3,043E-01 
3.043E-01 
3.041E-fl] 
3,043E-01 


Lb«/iec,  PRESSURE  IN 
f  HhP 

1.70PE-0?  1,501E- 

1.710E-02  1,501F- 

1.719E-02  1,504E- 

1.719E-02  l.SOOE- 

1.710E-02  1,498E- 

1.715E-fl2  I.SOIE- 

1.712E-fl2  1,495E- 

i.713E-02  1.496E- 

1.71lE-fl2  1.496E- 

1.71BE-02  1,:04E- 


Psia,  TEMPERATURES 
Hho 

■03  8,351E-03 
•03  8.351E-03 
■03  8.382E-03 
03  B.356E-03 
■03  8,375E-03 
03  8.355E-03 
03  8.342E-D3 
03  B.3B9E-03 
03  8,331E-03 
03  B.384E-03 


IN  R 
Phf 

19,94E4D1 

19.95E401 

19.98E401 

19,93E401 

19,90E+01 

19,95E401 

19,86E4D1 

19,88E401 

19.88E401 

19.98E401 


PhD 

91.20E400 
91.20E4CO 
91,546400 
91. 266400 
91.47E400 
91.25E400 
91,11E400 
9l.62Et00 
90,996400 
91,566400 


Count  Ta 

1  45,826401 

2  45,83E401 

3  45,836401 

4  45,83E401 

5  45.82''401 

6  45.B2E401 

7  45,826401 

8  45.82E401 

9  45,82E40I 

10  45,826401 


Tcin 

46.68E401 

46,696401 

46,696401 

46,70E401 

46,7IE4flI 

46,70E401 

46,69640! 

46,706+01 

46.706+01 

46,71E+01 


Texl 

15,506+02 
15.51E+02 
15,526+02 
15.526402 
15.516+02 
15,51E402 
I5,50E+02 
15  51Et02 
15,506+02 
15.50E+02 


Tex2 

14.026+02 

14.03E+O2 

14,03E+02 

14.036+02 

14.026+02 

14.026+02 

14.026+02 

14.036+02 

14.026+02 

14,016+02 


Tauqup 

90,186+01 

90,236+01 

90,256+01 

90,276+01 

90,296+01 

90,326+01 

90,346+01 

90,376+01 

90.396+01 

90.426+01 


Tauodl 

88,526+01 

88,54E+01 

88,576+01 

88,59E+01 

86,626+01 

e'’,65E+01 

88,676+01 

88.706+01 

88,736+01 

88,76E+81 


TaDQd2 
88.906+01 
8B.93E+0] 
88,956+01 
88,98E+D1 
89,016+01 
89,046+0  1 
89,076+01 
89,196+01 
89,116+01 
89,126+01 
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»*»*  PRE-RUN  DATA,  USING  AIR  ONLY  »»*» 

DATE;  3-tfi-n 
Pa=  580,8  Psia 
Hair=1.9013  (Lbn/sec) 

PPair  ,6157  (Ibn/sec) 

RUN  -  Air 


-  Air  heater  used 
No  catalvst 


«***  PRE-RUN  INPUT  »»«» 


Te5tno= 

BEHRENS 

Date= 

3-I1-90 

Fi)elid= 

Heaterfoei= 

HYDROGEN 

DairchoVe= 

.420 

Dhtfuchoke= 

.040 

Dhtcxchoke^ 

,0700 

Cdair= 

.970 

Cdh= 

.970 

Cdhor 

Gannahf= 

1.405 

Rhf= 

766,50 

K«bf= 

.1405 

Kfiho- 

,5589 

Knair= 

.5320 

Maird^ 

1,900 

Tcind= 

0.0 

Mhfd= 

.0028 

Mhod= 

.0228 

Hfueld= 

Xpa= 

33262,8800 

0pa0= 

.0817 

Kpc:: 

,3300 

844,8140 


»»*»  DATA  EXTRACTED  «««» 


Pc 

90,11E+P(! 
90,00E+00 
90,28E+0  0 
90,23E+00 
90,26E+0  0 
90,39E+0fl 
90,nE+00 
89,50E+00 
90,19E+00 
90,27E+00 


FUEL  FLOW  RATE  IN 
Count  Hair 

1  l,959E+00 

2  l,95i)E+(!fl 

3  l,953E+00 

4  1.955E+flfl 

5  1.951E+00 

6  1.949E4|)a 

7  l,952E+00 

8  l,953E+00 

9  l,95FE+00 

10  l,949£+00 

Count  Ta 

1  46,56Ef01 

2  46,57E+01 

3  46,57E+01 

4  46,57E+01 

5  46,57E+01 

6  46,57E+01 

7  46,57E+0I 

8  46,57E+01 

9  46,57E+0I 

10  46,57E+01 


Gal/nin,  GAS 
Hbpair 
6,343E-01 
6,314E-(I1 
6,323£-01 
6,331£-01 
6.316E-01 
6,3!0E-fll 
6,319E-01 
6,324E-0I 
6,341E-01 
6,3!3E-fll 

Tcin 

76.74Et01 
76,76E+01 
76,80E+01 
76.8(IE  +  01 
76,81E+0I 
76,76E+01 
76,76E+01 
76,77E+01 
76.77E+01 
76,80E+01 


FLOW  RATES  IN 
Mfuel 
1,650E-01 
t.65lE-Qt 
1,650E-01 
1,651£-01 
1,649E-01 
1,650E-01 
1,651E-01 
1.649E-01 
1,649E-01 
1,647£-01 

Te»l 

13,75EfD2 

13,75E+02 

13,74E+Q2 

13,74E^02 

13,71E+02 

13,74E+02 

13,72E+02 

13,72E+02 

13.71E+02 

13,71E+02 


Lb«/sec,  PRESSURE  IN  Psia,  TEMPERATURES 


9,(i96E-03 

9.t44E-03 

9,127E-D3 

?,118E-fl3 

9,133E-03 

9,144E-03 

9,I35E-03 

9.121E-03 

9.092E-03 

9.122E-03 

Tex2 

12  44E+02 
12.44E+02 
12,45E+02 
12,45E+02 
l2,44E+02 
12,44E+02 
12,44E+02 
12,44E+02 
12.44E+02 
12.44E«02 


Mhf 

3,821E 

3,823E 

3,823E 

3,824E 

3.824E 

3,B25E' 

3,82PE 

3,826E- 

3,B26E' 

3,825E- 


Mho 

Phf 

Pho 

1,941E-P2 

50.44E+01 

21,P6E*01 

t,940t-02 

SO.ATE^Ol 

2\,05Em 

1,937E-P2 

50,47E+P1 

21.D2E+01 

1.936E-P2 

50.48E+01 

2I,01E+01 

l,933E-02 

50  47E+01 

28,98EtDl 

1,928E-P2 

50,49E+01 

20.93E+01 

1,925E-P2 

50,4?E+01 

2D,e9E+01 

l,924E-02 

50,50E+01 

20.89E+01 

1.925E-P2 

50,50E+D1 

20,88E+01 

1.925E-02 

50,48Et01 

20,89E+01 

Taugup 

B9,34E40) 

B9,B9E+01 

89,89E+01 

89,91E*01 

89,92E+01 

89,93E+P1 

89,93E+01 

89,93E+01 

89,94E+01 

89.94E+01 


Taugdl 

88.46E+01 

88,49E+il 

B8.50E+01 

B8,51E+P1 

BB.52E+111 

8B,53^P1 

88,53E+01 

88,j5E+01 

BB,56E+01 

88,56E+01 


Taoqd2 

B8,75E+01 

8B,76E+01 

89,77E+01 

88,7BE+P1 

il’,^BPE+il 

88,81E+01 

88,83£tPl 

iim\ 
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««««  PRE-RUN  DATA,  USING  AIR  ONLY  «»*» 


DATE:  3-lt-90 
Pa=  577,2  Psia 
hair=l,8933  (LbK/sec) 
BPair  ,6131  (Ibw/sec) 


RUN  ^3  -  No  air  heater 

No  catalyst  used 


**«»  PRE-RUN  INPUT  «»» 


Te5lno= 

BEHRENS 

Da1e= 

Fuelid= 

Hesterfuel= 

HYDROGEN 

Dairchoke^ 

,420 

Dhtfuchoke= 

.040 

Dhtcxchoke= 

,0700 

Cdair= 

.970 

Cdh= 

,970 

Cdho= 

GafiNahf= 

1.405 

Rlif = 

766,50 

l(«hf= 

,1405 

X«ho= 

.5589 

XNair= 

.5320 

Maird= 

1,900 

Tcind= 

O.D 

Mhfd= 

.0028 

Mhod= 

.0228 

Mfueld= 

l(pa= 

33262.8800 

Vpa0= 

.0017 

i(pc= 

»»*»  DATA  EXTRACTED  «»»* 


eoo7 


,970 


,3300 

84VB140 


Pc 

90,53E+00 


90.53E+D0 


YU  . 

90,86E 


90,56E+00 

B9,92E+00 


90,56E+00 


90.77E+00 


90,49E+00 

90,56E+P0 


90,47E+00 


FUEL  FLOW  RATE  IN 

Gal/Min,  GAS 

FLOy  RATES 

Count 

Hair 

Mbpair 

Mfuel 

1 

1,8252+00 

5,911E-01 

3,2232-01 

0 

(. 

1.826E+00 

5,9142-01 

3,2422-01 

3 

1,8282+00 

5,9192-01 

3,2502-01 

4 

1,8252+00 

5.909E-01 

3.251E-01 

5 

1.8272+00 

5,9172-01 

3,2512-01 

6 

1.8332+00 

5,9372-01 

3.2532-01 

7 

1,8252+00 

5,9082-01 

3.2512-01 

8 

1.8252+00 

5.9092-01 

3,2502-01 

9 

1.8252+00 

5,9112-01 

3,2502-01 

10 

1,8252+00 

5,909E-01 

3,250E-01 

IN  Lbw/sec,  PRESSURE  IN  Puia, 
f  Hhf 

1,9071-02  1,DD0E-D4 

l,917E-02  l,000E-04 

1,920E-02  l,000E-04 

l,924E-02  l.OOOE-04 

l,921E-02  1,OOOE-04 

l,916E-02  l,D00E-04 

l,924E-02  1,000E-D4 

l,924E-02  l,000E-04 

l,923E-02  l,000E-04 

l,924E-02  l,000E-04 


TEMPERATURES  IN  R 
Mho  Phf 

1,0002-04  lD,0DE-05 

l.OOOE-04  10,00E-05 

l,00PE-04  lO.OOE-05 

1,OOOE-04  lO.OOE-05 

l,D00E-04  !0,OBE-05 

l,O00E-04  10,OOE-O5 

1,OOOE-04  lO.OOE-05 

1,OOOE-04  10,OOE-05 

l.OOOE-04  lO.OOE-CS 

1,0002-04  10,002-05 


Pho 

lC.OPE-05 
10.002-05 
10,OOE-05 
10,002-05 
lll,00E-05 
lO.OOE-05 
1  0,002-05 
lO.OOE-05 
1  0,0  02-05 
10,002-05 


Count 

Ta 

Tcin 

Teil 

1 

46.892+01 

47.83E+01 

16,942+02 

0 

L 

46.902+01 

47,832+01 

16.942+02 

3 

46,902+01 

47.842+01 

16,932+02 

4 

46,902+01 

47,842+01 

16.932+02 

5 

46,902+01 

47,852+01 

16.92E+02 

6 

46,912+01 

47,862+01 

16.922+02 

7 

46,922+81 

47,862+01 

16.922+02 

8 

46.922+01 

47,66E+01 

16.932+02 

9 

46.922+01 

47,862+01 

16.932+02 

10 

46,922+01 

47,872+01 

16,952+02 

Tex2  Taugup 
15.35E+02  99,58E+01 
15,35E+02  99,642+01 
15.352+02  99,68E+01 
15.34E+02  99.702+01 
15,35E+02  9?,71E+01 
15.35EtD2  99,722+01 
15.352+02  99.732+01 
15.352+02  99,732+01 
15.352+02  99,732+01 
15.362+02  99.732+01 


Taugdl 

96,782+01 

96,212+01 

95,732+01 

96.272+01 

96,792+01 

96,682+01 

95,692+01 

95,782+01 

96,382+01 

95,842+01 


Tarod2 
98.252+31 
98,282+01 
98,312+01 
98,342+0  1 
98.352+01 
96,382+0  1 
98,382+01 
98,392+0  1 
98.392+01 
98.402+0  1 
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»**«  PRE-RUN  DATA,  USING  AIR  ONLY  »««* 


DATE;  3-11-90 
Pa=  571,9  Psia 
Mair=i,8876  (Lbn/sec) 
Bhir  ,6112  (Ibn/sec) 


RUN  #4  -  No  air  heater  used 

Perlite  catalyst  installed 


»»»»  PRE-RUN  INPUT  »«»» 


Testno= 

Heaterfuel= 

BEHRENS 

HYDROGEN 

Date=: 

3-19-90 

Fuelid= 

Dairchoke= 

Dhtoichol(e= 

,420 

,0700 

Dhtfuchoke= 

,040 

Cdair= 

Ga«nahf= 

Rhf= 

Knhf=: 

,970 

1.405 

766,50 

,1405 

Cdh= 

.970 

Cdho= 

Ktiho- 

,5589 

Knair= 

.5320 

Mair(i= 

1.900 

Tcind= 

0,0 

Hhfd= 

,0028 

Mhod= 

.0228 

Mfueld= 

Kpa= 

33262,8800 

Vpa0= 

.0017 

i(pc= 

»»*»  DATA  EXTRACTED  **»» 


80  07 


,970 


,3300 

844,8140 


Pc 

94,21E+P0 

93,74E+00 

93,45E+00 


93,54E+00 

92,B5E+0C 


93,74E+00 

94,18E+00 


93,83E+30 

94,55E+0D 

93,56E+00 


FUEL  FLOW  RATE  IN  Gal/fiin,  GAS  FLOW  RATES  IN 
Count  Main  Mbpair  Mfuel 

1  l,912E+00  6,191E-01  3,178E-0l 

2  l,912E+00  6,192E-01  3,196E-0l 

3  l,916E+00  6,205E-01  3,201E-01 

4  l,912E+00  6,192E-01  3.206E-01 

5  l,910E+00  6,185E-01  3,207E-01 

b  l,912E+00  6,190E-Oi  3,210E-01 

7  l,916E+00  6.205E-01  3,211E-01 

8  l,713E+00  6,195E-01  3,210E-01 

9  l,913Et00  6,195E-01  3,210E-01 

10  1,918E+00  6.211E-D1  3,211E-01 


Lb«/5ec,  PRESSURE  IN  Psia,  TEMPERATUREE  IN  R 
f  Hhf  Mho  Phf 

l,795E-02  l,0D0E-04  l,0PDE-04  10,DOE-05 

1.805E-fl2  l,000E-04  l,000E-04  lO.OOE-05 

l,804E-02  1,OOOE-04  l,D00E-04  10,OrE-05 

l,811E-02  l,000E-04  l,000E-04  10,OOE-05 

l,8l5E-02  l,000E-04  1,OOOE-04  10,OOE-05 

1.814E-02  l,O00E-O4  l,O00E-04  l0,00E-05 

l,810E-02  l,000E-O4  1,OOOE-04  10,00E-05 

l,812E-02  l,000E-04  1,OOOE-04  10,OOE-05 

l,8l2E-02  l,0D0E-04  l,00PE-04  10.00E-t5 

l,80BE-fl2  1,OOOE-04  l,O00E-l)4  l0,D0E-05 


Pho 

10,OOE-O 

lO.OCt-0 

10,OOE-0 

lO.OOE-0 

16.00E-0 

lO.OOE-0 

10,O0E-0 

lO.OOE-O 

16.0PE-C 

10,OOE-0 


Count 

Ta 

Tcin 

Texl 

Tex2 

Taugup 

9B.13E+01 

Taugdl 

Tarad2 

1 

46,37E+01 

47.30E+C1 

15,53E+02 

14,14E+02 

72,85E+01 

71,3eEi01 

2 

46.38E+01 

47,30E+01 

15,53E+P2 

14,14E+02 

98.29E+01 

72,456+01 

71,39E+01 

3 

46,3?E+01 

47.30E+01 

15,55E+02 

14.14Et02 

98,37E+01 

72,96E+01 

71.42E+01 

4 

46,39E+01 

47,29E+01 

15.55Et02 

14.14E+02 

98.46E+01 

72,77E+01 

71.43E+P1 

5 

46,3?E+01 

47.29E+01 

15,53Et02 

14.13E+02 

98.53E+01 

72,85E+01 

71.43E+01 

6 

46,39Et01 

47.29E+01 

15.54Et02 

14.I4E+P2 

9B.59E+01 

72.21E+I1 

71,43E+01 

7 

46,38E+01 

47,29E+01 

15.55Et02 

14.15E+02 

98,67Et01 

71,B8E+01 

71.44E+01 

8 

46,37E+01 

47.29E+01 

15,55E+02 

14,16E+02 

98.76E+01 

71,76E+81 

71.45E+01 

9 

46,37E+01 

47,29E+01 

15.55E+02 

14.16E+02 

98.85E+01 

72.82E+01 

71,49E+01 

10 

46,36E+01 

47,28E+01 

15.55E+02 

14.17E+02 

98,93E+01 

73,24E+81 

71,50E+01 
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